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ABSTRACT 

The film-c:ondsnsatioa character isuics of a smooth tube 
and six externally-finned tubes having fins 1 tnn high and 1 
mm thick, and pitches of 1.5, 2.0, 2.5, 3.0, 5.0, and 10.0 

mm, were experimentally tested. 

A smooth copper tube with an active length cf 133.5 mm, 
an outside diameter of 19.05 mm, and an inside diameter of 
12.7 mm was first tested to correlate the inside heat- 
transfer coefficient using the Sieder-Tate equation. The 
leading coefficient for this equation was found to be 0.034 
± 0.001, and was used to derive the excernal condensing 
coefficient for all of the tubes by subtracting the inside 
and wall resistances from the measured overall resistance. 
The condensing coefficient was measured, both ar atmospheric 
pressure and vacuum (84 mm Hg) , wirh the heat flux as a 
variable. 

Condensation data taken for the smooth tube were 
compared with data in the literature to check the reli- 
ability of the apparatus and the data-reduction procedures. 
The data for the finned tubes showed an optimum pitch of 2.5 
mm. 
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NOMENCLATURE 



c p Specific heat, JcJ/kg-K 

Dj Inner diameter, m 

Do Outer diameter, m 

F Parameter, g(^hf^/u2kZf 

g Acceleration due to gravity, m/s^ 

Gr^ Grashof number of air, gO(T*^-T^) 

h Heat-transfer coefficient, ‘rf/m^-K 

hfj Latent heat of vaporization, kJ/kg 

k Thermal conductivity, w/m-K 

k^ Thermal conductivity of air, W/m-K 

kg. Thermal conductivity of glass, W/m-K 
k^ Thermal conductivity of liquid, W/m-K 

L Length, m 

Lg Length of the toiler, m 

Nu Nusselt number, hd© /k 

p Pressure, MPa 

Pr Prandtl number, Cp/k 

Q Heat transfer rate, W 

g Heat flux, H/m^ 

R Thermal resistance, K/W 

Ryv Wall thermal resistance (D^^ln (D^, /D,- ) /2k^) 



K/W 
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r,-^3 Inner radius of the boiler, m 

r<3 3 Outer radius of the boiler, m 

Ra Rayleigh number, GrPr 

Re Reynolds number, u^Do 

Re Two-phase Reynolds number, 

T Temperature, ®C 

T^ Film temperature, °C 

T;^g Inner wall temperature of the boiler, ° C 

Tj,.g Outer wall temperature of the boiler, °Z 

T« Ambient temperature, ®C 

T Average temperature, °c 

u Velocity, m/s 

u^ Vapor velocity, m/s 

Oo Overall heat- transf a r coefficient, W/m^-K 

X Parameter, D-^^i/ 3 /Zk^ 

Y Parameter, 

Z Sieder-Tats parameter, Re^- ®Pr^ (/'/^^) o • 

Q Expansion coefficient of air, 

^ Dynamic viscosity, N-s/m2 

Kinematic viscosity, m 2 /s 
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I. INTRODOCTIDN 



A. BACKGROONO 

As combat systems aboard naval surface ships continue to 
grow in importance, size, and weight, the feasibility of 
their installation aboard lightweight vessels can only be 
realized by reducing other major shipboard weight require- 
ments. A significant achievement in reducing main propulsion 
weight was reached with the advent of gas-turbine propulsion 
plants, of which the majority of our new combatants will be 
powered. However, for the remaining capital ships and subma- 
rines to be constructed using either conventional or nuclear 
steam systems, power plant weight reduction can only be 
realized by increasing the effectiveness of the individual 
plant components and thereby reducing their size. Plans call 
for future installation of Rankine-o ycle waste-heat recovery 
systems even on those ships with gas-turbine propulsion, so 
significant weight reduction could also be realized there 
through the use of more effective systems. 

Due to the nature of the condensing process, the 
greatest thermal resistance to condenser-tube heat transfer 
occurs on the tube side, and a thorough review of tube-side 
enhancement is presented by Bergles and Jensen [Ref. 1 ]. 
Internal enhancement of the tube without regard to the 
shell-side resistance problem, however, would be a wasted 
effort if improvements were made to the point where the 
external (shell-side) resistance became the controlling 
factor in the transfer of heat. Recent investigations into 
the shell-side condensing process and techniques to improve 
this process are thoroughly reviewed by Marto [Ref. 2]. 
Outside enhancement techniques include the use of low 
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integral fins, roped tubes, fluted tubes, and applied 
coatings to promore dropwise condsnsa rion. 

Ongoing research is underway at the Naval Postgraduate 
School to further study the shell-side enhancement of marine 
condensers. An endurance apparatus is in operation to 
examine the dropwise-promoring effectiveness of various 
metallic and polymer coatings over prolonged periods of 
time. Reilly [ Bef . 3] conducted a study of improvements 
using various spirally fluted tubes under filmwise condensa- 
tion, and the effects of tube bundle inundation have also 
been studied by Kanakis [Ref. 4]. 

A test apparatus has been constructed by Krohn [Ref. 5] 
to systematically study steam condensation on a single hori- 
zontal tube. The apparatus is modeled after a similar design 
constructed by Rose [Bef. 6], but can also operate under the 
high-vacuum conditions found in marine condensers and 
without the presence of noncondensable gases. This test 
apparatus was instrumented and tested by Graber [Bef. 7], 
who determined the Sisder-Tate coefficient for the test tube 
length to be 0.029. 

In experiments using low integral fins at Queen Nary 

College of London, Rose, et. al. , [Ref. 6] obtained an 

optimum pitch of 2.0 mm (0.08 in.) for a gnven fin size of 1 

mm (0.04 in) high and 0.5 mm (0.02 in) wide for tests 

conducted at one atmosphere. with this pitch, he reported 
vapor-side heat-transfer enhancements of 400^ and 300% for 
heat fluxes of 0.3 MW/mz (9.51x10^ Btu/hr-ftz) and 0.8 MW/mZ 
(2.54x105 3tu/hr-ftz), respectively, and for a vapor 
velocity of 0.7 m/s (2.3 ft/sec). Further testing is needed 
to compare these results with similar fin geometries under 
vacuum conditions, while varying fin dimensions (height and 
width) in addition tc the fin pitch. 
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B. OBJECTIVE 



The main objectives of this research effort were, there 
fore, to: 

1) ensure a vacuum -tight apparatus so that data could be 
taken at both atmospheric and vacuum conditions with no 
detrimental effects due to the presence of nonconden- 
sable gases; 

2 ) take data for an instrumented smooth tubs to verify 
the Sieder-rate coefficient obtained by Graber; 

3) take data for a smooth tube to check the reliabilty 
of the apparatus and the dat a-redu ctio n procedures used; 
and 

4) take data for six externally-enhanced tubes of 
various fin pitches to obtain the relative optimum pitch 
for a fixed fin geometry. 
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II. DES CRIP TION OP APPA RAT US 



A. STSTEH OVERVIEW 

The apparatus used for this research was essentially the 
same used in references 5 and 1 , with several noted modifi- 
cations. A schematic sketch of the system is shown in 
Figure 2.1. Steam was generated in a 30U.8-mm (12-in) diam- 
eter Pyrex glass boiler by ten 4003-watt, 480-volt Watlow 
immersion heaters. Passing througa a 304.8-mm (12-in) to 

152.4-mm (6-in) reducing section, the steam travelled upward 
through a Pyrex section 2.44 m (8.3 ft) in length, around a 
180-degree bead, and back down a straightening section 1.52 
m (5.0 ft) in length before entering the stainless-steel 
test section. The condenser tube to be tested was mounted 
horizontally in the test section behind a viewport to permit 
visual observation of the condensing process. 

Steam that did not condense on the test tube passed into 
a stainless steel auxiliary condenser, and all condensate 
was returned via gravity to the boiler. The auxiliary 
condenser was constructed of two 9.5-mm (3/8-in) water- 
cooled copper lines helically coiled to a height of 457 mm 
(18 in) . 

Cooling water for the test tube was provided by two 
centrifugal pumps connected in series. The water could be 
throttled from zero flow to 0.69 1/s (11 gpm). The maximum 

water velocity which could be obtained through the tube was 
5.48 m/s (18 ft/sec) . A continuous supply of tap water was 
used for cooling the auxiliary condenser. Throttling the 
flow of tap water through the condenser was the means used 
to vary the internal pressure of the test apparatus. The 
water flow through both the test tube and the auxiliary 
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Figure 2.1 Scliematic of the test apparatus 




condenser was regulated by 19.1 aiin (3/4 in) dianeter needle 
valves and measured by rotameters with full-scale ranges of 
0.69 1/s (1 1 gpm) . 

An air ejector provided for removal of non condensable 
gases from the auxiliary condenser through a 12.7-mm 
(1/2-in) line. The source for the air ejector was 1.1-MPa 
(160-psig) house-air supply. 

B. SYSTEM INSTHOMENTATION 

The input voltage through the heaters was varied through 
a panel-mounted potentiometer. 440-VAC line voltage was 
reduced by a factor of 100 when fed into a differencial 
input precession voltage attenuator. The stepped-down 
voltage passed through a True-Root-Mean-Square convercer 
stage on which the integrated period was reduced to about 1 
ms. The output of the TRMS converter was then buffed and 
compared to a reference voltage from the potentiometer. The 
comparator output was fed to the control input of a Halaar 
silicon-controlled rectifier power supply which applied the 
actual volcage to the heaters. The IRMS converter output was 
also paralleled to a filter and than input zo the data 
acquisition system. This input was proportional to the p'ower 
supply output. A diagram of the system is shown in Figure 
2 . 2 . 

The internal pressure of the system was measured manu- 
ally by a U-tube, mercury -in-g lass manometer graduated in 
millimeters. Unavoidably, steam could condense in the manom- 
eter. Therefore, the varying height of the water column in 
the manometer needed to be accounted for when measuring the 
system pressure. 

Temperatures throughout the system were measured by 
copper-constant an thermocouples: six for the wall of a 

specially-constructed test tube, two for the steam, and one 
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Figure 2.2 Line diagrau of the power supply. 



each for the cooling-water inlet, condensate return, and 
ambient. The calibration procedure for these rh ermoccuples 
is described in Appendix A. The temperature rise through the 
tesn tube was measured by a Hewlett-Packard (HP) 2804A 
quartz thermometer. 

All temperature measurements were fed directly into the 
data-acquisition system as described below. 

C. DATA ACQOISITION 

An HP 3497A Data Acquisit ion/Control Unit was used to 
monitor system temperatures. This was interfaced with an HP 
9826A computer which served as a controlling unit through an 
interactive da ta-reduction program and user keyboard 
prompts. Saw data gathered by the data-acquisition system 
were stored on computer disks for later reduction and evalu- 
ation . 

D. SYSTEM HODIFICATIONS 

1 . B oile r 

The fiberglass insulation was removed from the 
boiler to allow the operator to acre easily monitor the 
water level. Although a closed-system design was used, it 
was still possible for steam to escape via the air ejector 
or through the relief valve [Fig. 2.1]. Calculations showed 
that the additional heat loss due to the removal of the 
insulation was minimal [App. 3], and the author felt this 
loss was much more acceptable than risking damage tc the 
immersion heaters through a low-water casualty in the 
boiler. 
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2 • Con den sa_^ ? ipi nq 

As originally designed, draining rhe system required 
breaking down the condensate piping. While this was not a 
daily cccurrance, the procedure was inconvenient and there 
existed also the possibility of losing the vacuum integrity 
of the system each time it was done. 

To avoid these problems, the existing fill-line 
valves were rearranged as shown in Figure 2.1. This 
arrangement also added two additional features to the 
system: 

1) the fill/drain valve could be opened during operation 
to drain any heavy particulate matter from the system - 
similar to the "bottom blow" procedure used on Naval 
boilers; and 

2) after extended periods of inactivity while opened to 
the atmosphere, the entire system could be given a thor- 
ough steam-cleaning by following the procedures outlined 
in Appendix C. 

2* Vent Valve 

The modification of the condensate return piping 
necessitated the addition of a vent valve for use when 
filling or draining the system. A 4.3-mm (0.17-in) needle 
valve was installed on the 101.6-mm (4-in) flange of the 
test section. This valve would also serve as the tap for the 
proposed sampling of nonco ndensabie gas concentrations in 
the system. 

4 . Man ome ter Li ne 

The original system design used a 6.4-mra (1/4-in) 
stainless steel tube angled down to the mercury manometer. 
During this thesis, the manometer was raised to eye level to 
facilitate easier and quicker reading. Replacement of the 
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srainlsss-steal line by a 12.7-ina (1/2-in) copper tuba 

reduced the possibilty of error caused by water slugs 
building up in the smaller diameter tube. The more workable 
copper was chosen over stainless steel to reduce the stiff- 
ness of the connecting line. This was necessary to 
eliminate leakage in this part of the apparatus as explained 
in Section II. E. 

5 • P ress ure Tra n sducs r 

As an alternative to the manometer, a Celesco 
strain-gage pressure transducer was installed on the test 
section flange next to the vent valve. The calibration line 
for the transducer is shown in Figure 2.3. The author felt, 
however, that the reliabili/ty of this measurement would not 
be high enough until a second, more accurate transducer was 
installed. Once incorporated into the system, though, these 
transducers would provide automatic input to the data acgui- 
sition system, eliminating the requirement to manually enter 
the manometer reading into the data- r educt ion program. 

6 • Relief V alve 

As originally designed, a 5895-Pa gage pressure 
(1.0-psig) relief valve was installed beneath a 1.0-m 
(39.4-in) length of 12.7-mm (1/2-in) stainless-steel piping. 
Steam which condensed and became trapped in this piping 
would open the valve at a water-column height of only 0.70 m 
(27.6 in). Once opened, a back-pressure of 0.14 HPa gage 
(20 psig) was required to reseat the valve - something unob- 
tainable even with an absolute vacuum on the inlet side of 
the valve. To avoid this problem, the valve was raised in 
the line to a point only 7 6 mm (3.0 in) below the outlet 
from the auxiliary condenser section and was replaced by 
another 6895 Pa gage pressure (1.0-psig) relief valve which 
reseated at only 0.04-MPa gage (6-psig) back- pres sure. 
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Figure 2.3 Calibration of the pressure transducer. 
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St g am- te m per atur a Probe 

The steam temperature prooe was located directly 
above the test tube and shared tha same inlet to the test 
section as the manometer line. When drawing a high vacuum, 

this arrangement allowed water which had collected in the 
manometer to be drawn back into tha test section where it 
would flow down the probe and onto the tube. To prevent this 
water and any contaminants picked up in the manometer from 
being deposited on the tube, the probe was bent sc that it 
was offset from the center of the test tube. 

8- Co olin g Water System 

A second centrifugal pump was added in series to the 
one already installed. This pump boosred the maximum cooling 
water velocity through the test tube from 3.96 m/s (13 
fr/sec) to 5.48 m/s (18 ft/sec) . 

The 10.55 kW (3 Ton) air-conditioning unir used with 
the cooling water system would energize at a water tempera- 
ture of 17 ®C (62.6 ^ F) and secure when the temperature was 
reduced to 13.4 °C (56.1 ‘’F). Therefore, for a given steam 
temperature of 50 ®C (122 °F) around the tube, the log- 
mean-temperature difference would vary by as much as 11%. 
But the measured temperature rise of the water through the 
tube showed very little change. To avoid rhis rransient 
problem, the air-conditioning unit was not used, and instead 
fresh tap water was continuously fed to the sump while an 
equal amount of water was being drained, maintaining a 
constant sump level. This method provided a constant- 
temperature supply of cooling water to the inlet side of the 
test tube. 

It should be pointed out that, while the duty cycle 
of the air-conditioning unit was a function of the ther- 
mostat used, a more sensitive thermostat would require the 
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US 9 of a more expensive cooling device than a commercial air 
conditioner. An alrernative solution would be to install a 
much larger sump in the system. 

9 . Thermo pile 

Since the temperature rise of the cooling water 
through the tube was a critical measurement in the experi- 
ment, a 10-junction thermopile was added to measure this 
temperature rise in addition to the quartz thermometer. As 
will be explained in Section V.A.3, however, problems arose 
with the thermopile, and it couli not be accurately used 
during this thesis. 

E. VaCOOH INTEGRITT 

One of the objectives of this thesis was to ensure a 
vacuum-right test apparatus to elimlnare rhe presence of any' 
noncondensable gases and their detrimental effects on the 
condensing process. A standard of no more than a 5.0-mm 
mercury (0. 10-psi) loss over a 24-hr period was considered 
to be an acceptable tolerance, but obtaining a leak rate 
within this tolerance proved to be the most time-consuming 
effort during the research. 

Due to the construction of the apparatus and nature of 
the experiment, most leak -detection methods could not be 
used. A sealing substance could not be used without risking 
contamination of the interior of the apparatus which might 
prohibit filmwise condensation on taa test tube. 

Initially the system was pressurized and the standard 
soap- solution test was used to locate leaks. Once pressur- 

ized, a liquid-soap solution was applied to each external 
joint or fitting where a leak could be present. The higher 
pressure air inside the apparatus would escape through any 
leaks and produce bubbles on the applied soap film. 
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Howsver, the maximum pressure which the Pyrex glass members 
could tolerate was only 0.074 MPa gage (10.7 psig) , and for 
safety reasons the author chose not to pressurize the system 
to more than 0.034 MPa gage (5.0 psig) - a pressure differ- 
ence between the atmosphere and the apparatus of only 0.040 
MPa (5.8 psi). The numerous external valves and fittings 
prohibited the use of an evacuated hood to achieve a greater 
pressure difference. 

A similar test could not be used to locate any vacuum 
leaks which were not present when the system was under a 
positive pressure, as the apparatus was not large enough to 
permit the application and observation of a soap solution on 
the interior. 

In another attempt to locate the leaks, a National 
Research Corporation (NRC) 101.6-mm (4-in) vacuum pumping 

system was connected to the apparatus. This pumping system 
included a Welch model 1376M mechanical pump and a model 
NHS-4 diffusion pump. An NRC model 521 thermocouple gage was 
also connected to the test apparatus and the entire appa- 
ratus was evacuated to 0.21 torr (4.1x10~3 psia) . Acetone 
was sprayed around all flanges, fittings, and joints. A leak 
around any of these should have produced a rapid rise in the 
thermocouple reading, but this method also proved ineffec- 
tive, probably due to the large size of the apparatus 
resulting in too great a mean free path for the acetone 
molecules to travel from the leak to the thermocouple. 

The next alternative was to break the system apart into 
three main sections; the glass boiler and steam piping, the 
stainless-steel test section and auxiliary condenser, and 
the condensate return piping. The glass section was blank- 
flanged and evacuated to an absolute pressure of 0.033 torr 
(6.4x10 -♦ psia) . The rate-of-rise measurement for this 
section showed a loss of only 0.43 mm of mercury (0.01 psi) 
over a 30-hr period. It was, therefore, concluded that any 
leak in the assembled apparatus was not from this section. 
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Once removed, the test section and the auxiliary 
condenser were blank-flanged, pressurized to 0.10 MPa gage 
(15 psig) , and immersed in a large plsxiglas tank filled 
with water. This test easily revealed a number of small 
leaks abcut the inlet side of -he tesr tube and also around 
the plug which connected the condensate return piping to -he 
base of the auxiliary condenser. Replacing an 0-ring in the 
tube fitting and silver brazing the plug into place elimi- 
nated these leaks. 

The same immersion test revealed small leaks in the 
joints of the condensate return piping. These leaks were 
eliminated by replacing all stainless-steel ferrules in the 
Swagelck fittings with teflon ferrules. 

Once reassembled, considerable leakage was still indi- 
cated by a substantial overnighr rise in the manometer 
level. The author felt confident that this leak was not in 
the main assembly of the apparatus, but was instead in the 
manometer assembly itself since it could not be pressurized 
for testing. As mentioned in section II. C. 4, the stainless- 
steel line leading to the manometer was at this time 
replaced by a 12.7-mm (1/2-in) soft copper tube. This tube 
eliminated the need for two 90-degree elbows and three 
lengths of stainless steel tubing in the line - an assembly 
which proved too rigid to allow even the slightest misalign- 
ment into the manometer. 

Open completion cf the installation of this assembly, 
the system was evacuated to an absolute pressure of 92.5 mm 
Hg and ever a 24-hr period the mercury level rose to only 
94.0 mm. This leak rats was well within the acceptable 
tolerance. 
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F. TOBES TESIED 



1 • Ins tru m ent ed Tube 

An instrumented tube was fabricated from a thicic- 
walled ccpper tube with an inner diameter of 12.70 .nn (1/2 
in) and an outer diameter of 19.05 mm (3/4 in). The tube 
was cut inro three sections inno which six holes were 
drilled axially along the walls at equal spacings 60 apart. 
These passages were fitted with 0.094-mm (3/32-in) OD capil- 
laries [Fig. 2.4] which were silver-soldered into place, and 
the three sections of tube were rhan soldered back into one 
piece. Thermocouples were fitted into the capillary 
sections to measure an average wall temperature. 

By knowing the average wall temperature, the Nusselt 
number for the inside could be compu-ad. By computing the 
gradient of the Musselt number againsz the Sieder-Tate 
parameter, the inside coefficient could be obtained as -he 
inverse of the gradient. Figure 2.5 shows a photograph of 
the instrumented tube with the installed thermocouples. 
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Figure 2.4 Schematic of the instrumented tube construction 
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Figure 2. 5 Photograph of the instrumented tube. 
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2. Sm ooth Tube 

A smooth tube with no wall thermocouples was also 
tested to obtain the inside heat-transfer coefficient 
through the use of a modified Wilson Plot [Ref. 8]. 
Determination of the inside heat-transfer coefficient was 
critical to the experiment, as it was used to obtain the 
outside condensing coefficient for the smooth tube and all 
of the enhanced tubes. 

3 • Finne d Tub es 

To fulfill the main objective of this thesis, a 
series of six finned tubes was also tested [Fig- 2.6]. 
These tubes had the same overall dimensions as those above, 
but were enhanced with radial fins 1 mm (0.04 in| high and 1 
mm (0.04 in) thick. Each tube had a different fin pitch and 
was tested to determine a relative optimum pitch. Fin 
pitches tested were 1.5, 2.0, 2.5, 3.0, 5.0, and 10.0 mm. 

G. SYSTEH OPERATION 

The tube to be tested was cleaned in a warm solution of 
Sparkleen and then rinsed with tap water, which produced a 
contaminant-free, wetted surface. The tube was then 
installed in the test section, cars being taken not to touch 
or contaminate the condensing surface. 

The system was brought to operating pressure by 
following the procedures of Appendix D, and data collection 
began when steady-state conditions were achieved. 
Steady-state conditions were determined by observing the 
steam temperature measured by the respective thermocouples. 
When their output voltage on the HP 3497A reached a constant 
value with fluctuations of only one or two microvolts, it 
was assumed that steady-state conditions existed in the test 
apparatus. 
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Figure 2.6 Photograph of the finned tubes tested. 




Data sets were taken starting wirh a test-tube cooling 
water flow rate of 30% (which corresponded to a water 
velocity of 4.95 m/s or 16.2 ft/sec} , ranging downward in 
decrements of 10% through a minimum flow of 20% (1.16 ra/s, 

3.8 ft/sec), and then upward from 25% to 85% in increments 
of 10%. After adjusting the flow rate, the temperature rise 
through the tuba was monitored by observing the digital 
output of the quartz thermometer. When this rise became 
constant, another set of data could be taken. During data 
runs, a slight rise in pressure accompanied the decrease in 
the cooling water flow through the tube, a result of the 
reduced heat flux. Similarly, increasing the flow through 
the tube caused a slight decrease in the system pressure. 
This variance could be anticipated, and since data were to 
be collected at a constant pressure, it was easily compen- 
sated for by throttling the flow of coolrng water through 
the auxiliary condenser one or two percentage points on the 
rotameter. 

Something which could not be anticipated, however, were 
sudden fluctuations in the tapwater pressure to the auxil- 
iary condenser which caused pressure changes of several 
millimeters of mercury in the system. To avoid this 
problem, the flow through the auxiliary condenser had to be 
continuously monitored, unless data were being taken lata at 
night when there was no demand on the laboratory building's 
water supply. The test tube could be easily monitored 
through the viewport for confirmation of filmwise condi- 
tions. If there was any sizeable change to dropwise 
condensation on the tube, the data set was discarded and the 
procedure was repeated. 
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III. FILMWISE CONDENSATION 



A. THE DROPHISE PROBLEM 

It was essential during the course of this thesis to 
collect data under filmwise conditions. Numerous problems 
were encountered by Graber (Ref. 7] in avoiding the tran- 
sition zo dropwise condensation during operation, and his 
proposed solution wa s a vacuum-tight test apparatus. Since 
the tube surface would wet completely after installation, 
contaminants leaking into the system were possibly adhering 
to the surface and promoting the dropwise condensation. 

However, even after obtaining a vacuum-tight apparatus, 
the author was still unable to maintain good filmwise 
condensation for more than two hours on a smooth tube. While 
this was enough time to collect a complete set of data for 
this tube, filmwise condensation lasted seventeen minutes at 
the most on any of the finned tubes, the average time being 
less than ten minutes. This was predictable - the corners of 
the fin/surface interfaces provided a better trap for 
contaminants and were harder to clean - but unacceptable. 

The use of the steam-cleaning method described in 

Appendix C would thoroughly clean ~he tube so that complete 

filmwise condensation was r e- established, but dropwise 

condensation would again become prevalent within minutes. 

B. SOLOTIOH 

Having eliminated the possibilities of iastalling a 

dirty tube or contamination due to leakage, the only reason 
for the dropwise problem had to be coming from outgassing of 
the nylon holders for the test tube. The outgassing rate for 
nylon was found to be almost two orders of magnitude greater 
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than the rate for Teflon [Ref. 9], Teflon being the same 
order as stainless steel. By puiping down the apparatus 
while it was hot, nylon molecules were being outgassed into 
the test section and were immediately being deposited onto 
the surface of the cooler test tube, resulting in inadver- 
tent "sputtering” of the tube with a nylon coating and 
subsequent dropwise condensation. To eliminate the nylon 
interface with the interior of the test section, special 
stainless steel caps were manufactured to slip over the 
nylon holders. 

Teflon bushings were fitted within the caps to insulate 
them from the the tube. Teflon had both a Low thermal 
conductivity and a relatively low outgassing rate. Figure 
3.1 shows a detailed sketch of the installation. 

This configuration appeared to solve the problem of 
dropwise condensation. Although an actual endurance test was 
not conducted, rhe system was in operation intermittently 
for over fifteen hours with the smooth rube and over four 
hours with each finned tube with no breakup of the filmwise 
cond ensation. 

Once installed, this arrangement also eliminated the 
need for the tube-cleaning procedure recommended by Graber 
[Ref. 7], who felt that a strong cleaning solution of sodium 
hydroxide and ethancl was needed to decontaminar e the 
surface. Only a warm solution of Sparkleen was used 
throughout the data-ccllection stage of this thesis. 
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Figure 3.1 Cross-sectional view of the tube and holders. 
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VaPOfi VELOCITT LIMITATIONS 



To check the outside hea t- transf er coefficieat data with 
the Nusselt prediction, and also to obtain an accurate 
inside coefficient with the modified Wilson Plot, vapor 
velocities approaching zero were preferred. Due to the 
design of the system, however, the pressure drop to the 
auxiliary condenser required vapor flow past the test tube. 
This being the case, attempts wars made to minimize this 
flow velocity by cutting down the power to the boiler and 
throttling back on the cooling water supply to the 
condenser. As vapor velocity was decreased, however, drop- 
wise condensation again took place on the test tube. Under 
atmospheric pressure, this occurrence was at vapor veloci- 
ties of about 0.5 m/s (1.6 ft/sec), and under vacuum 
operations it occurred at about 0.9 m/s (3.0 ft/sec). 
Apparently there still existed a sizeable rate of outgassing 
within the test section, those gases collecting about the 
test tube and interfering with the filmwise condensation 
process. This suspicion was confirmed when an increase in 
the vapor velocity eliminated the formation of drops on the 
tube . 
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IV. DATA REDUCriON 



The data collected and stored on the computer disks 
reduced using the programs WILSON, SIEDER, and DRP. The 
programs were amenable to changes, which allowed the aurhor 
to analyze and compare results while varying parameters 
within the programs. Stepwise reduotion procedures for each 
program are listed below and the program listings are found 
in Appendix E. 
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7. B ESOLT S MS discussion 



Numerous data runs were made using the procedure 
described in Section II. G. Time constraints, however, 
limited the number cf repeat runs than could be made for 
this thesis. Primary concern focused on establishing a reli- 
able, repeatable Sieder-Tate coefficient. Dana were taken 
for all of the tubes at both atmospheric and vacuum (88mm 
Hg, 1.7 psia) conditions. Complete filmwise condensation was 
maintained for all data runs, and the mass concentration of 
noncondensable gases was held between 0% and during all 

testing. The negative value was indicative of slight super- 
heat in the system or an inaccurate manometer reading. 

A. INSIDE HEAT TRANSFER COEFFICIENT 
1 • Ins tru me nted Tube 3 esul ts 

Figure 5.1 shews the variation of the Nusselt number 
as a function of the Sieder-Tate parameter for the instru- 
mented tube run at atmospheric pressure (run S0A7) . This 
method yielded a Sieder-Tate coefficient of 0.035 on two 
separate data runs (SDA7 and SDA8) . The same method under 
vacuum conditions yielded a coefficient of 0.037, which was 
also repeated (runs SDA5 and SDA6) . The temperature distri- 
bution around the tube wall was symmetrical about the 
vertical plane passing downward through the centerline of 
the tube, and showed as much as a 15 C temperature drop 
from the top of the tube to the bottom. 

2. Smooth Tube R esult s 

Figure 5.2 shows the modified Wilson Plot for 
smooth-tube data collected at atmospheric pressure (run 
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Figure 5.1 Inside Nusselt nuiaber plot for the instrumented tube. 
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Figure 5.2 Hilson plot for the smooth tube. 
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0D12) . This method produced a Sisder-Tate coefficient of 
0.033 with an intercept of 1.16. A similar plot for vacuum 
conditions (run OD10) produced a coefficient of 0.036 and an 
intercept of 1.27. The Nusselt theory predicts an intercept 
of 1.53 but, as was noted in Section III.C., a finite vapor 
velocity past the tube was unavoidable in this research. 

Figure 5.3 shows a plot of the Fujii [ Hef . 10] 
correlation for the smooth-tube runs. This experimental 
correlation is to correct the Nusselt number for the effects 
of vapor shear on the test tube and is a plot of the equa- 
tion : 

NuRe-»/2 = 0.96FI/5 

The results obtained were slightly higher than those 
predicted using this correlation. 

3 • ere pa ncie s 

The original Sisder-Tate equation [Ref. 11] for 
fully developed turbulent flow in a tube with an L/D ratio 
greater than 60 has a leading coefficient of 0.027, sc a 

higher value for the shorter tube (L/D = 18) used for this 
research stands to reason. The data for the tabes appears 
at first inconsistent, with results of 0.033, 0.035, 0.036, 

and 0.037 - showing up to a 6% scatter from the mean. Both 
tubes, however, showed a larger coefficient when the test 
section pressure was reduced, which reduced the heat flux as 
well. Dropping the pressure from atmospheric to a high 
vacuum (88 mm Hg, 1.7 psia) reduced the heat flux by a 
factor of almost three [Fig. 5.4]. Stated another way, a 
low heat flux produced coefficients of 0.036 and 0.037, 
while a higher heat flux produced 0.033 and 0.035. 

The author doubts the reliability of the inside co- 
efficient obtained under vacuum conditions for two reasons: 
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Figure 5.3 Fujii correlation [ Bef. 10] for smooth-tube data. 
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Figure 5.4 Relative heat fluxes for atmospheric and vacaum runs. 



First, reducing the flow rate through the tube 
lowers the heat flux and increases the inside convective 
thermal resistance. Since the condensing coefficient is 
inversely proportional to the heat flux raised to the one- 
third power, the outside condensing resistance decreases. 
Because of this, the inside resistance becomes the dominant 
factor in the overall thermal resistance of the tube. Any 
error in measuring the inside heat transfer coefficient will 
greatly amplify the computed error in the condensing resis- 
tance. Thus, the condensing coefficient is highly sensitive 
to the accuracy of the inside coefficient, particularly at 
low water velocities (and corresponding heat fluxes) . 

Second, the test apparatus was subject to radio- 
frequency (RF) interference from the power supply. As the 
power was cut bade to lower the heat flux, the silicon 
controlled rectifier chopped the input signal proportionally 
and emitted the chopped portion as RF energy. Despite 

efforts to shield the cabling to the data- acquisition 

system, it still captured the RF signal with unacceptable 

results, and the thermopile was experiencing instantaneous 
fluctuations of several hundred microvolts. Increasing the 
power supply, however, decreased the errant signal until, at 
the maximum power input used for the high heat flux measure- 
ments, the thermopile showed a near-steady temperature 
difference that was within 0.1 ®C of the quartz thermometer. 

For these reasons, the 3ieder-Tate coefficients 

obtained from the low heat-flux runs were discarded. 
Averaging the coefficients found from both measurements for 
the high heat flux case yielded a coefficient of 0.034 ± 
0.00 1 . 
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FINNED TUBES 



Figure 5.5 shows a relative plot of the outside 

condensing coefficient for all six finned tubes compared 
with smooth tube results and the Nussalt line. Figure 5.6 
shows the same results corrected for area increases due to 
finning. These data were computed using the leading coeffi- 
cient on the inside of 0.03U obtained in the previous 
seer ion. 

All tubes tested show good enhancement over the smooth 
tube. Notice the increase of the smooth tube over rhe 
Nusselt line - a result of the inherent vapor shear in the 
test section. The increased scatrer in the lower heat-flux 
range is again a function of the increased significance of 
errors in measuring the inside thermal resistances. The 
plotting program used neglected ordinate values below zero 
or greater than IQs, so the scattering in this portion of 
the plot is actually worse than it appears. Figure 5.7 is a 
plot of the values obtained for the outside heat transfer 
coefficients of the six finned -ubes for a constant heat 
flux of 250,000 W/m 2 (79,000 Btu/hr-ft^) and a pressure of 
88 mm Hg (1.7 psia). This plot more clearly shows the 
optimum pitch of 2.5 mm. 

The heat-transfer characteristics for any tube will be 
enhanced by the addition of fins. In the case of purely 
convective heat transfer, the enhancement is a function of 
the increased surface area exposed to the fluid medium. 

Filmwise condensation, however, is dissimilar in that 
the build-up of a condensate film acts as an additional 
thermal resistance for the hea t-transf ar process. The objec- 
tive of tube enhancement, therefore, is to decrease this 
film thickness. The surface tension forces in the conden- 
sate tend to draw the liguid toward the fins, leaving the 
tube surface with a thinner film. The thinner film results 
in a higher heat -transfer coefficient. [Bef. 12]. 
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Figure 5.5 Co*parison of finned tubes with smooth-tube performance. 
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Figure 5.6 Outside condeasing coefficients corrected for area. 
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Figure 5.7 Variation of condensing coefficient with fin pitch 
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The optimum pitch found in this research increases nhe 
surface of the condenser tube 88%, but condensing heat- 
transfer coefficients for that pitch were enhanced by as 
much as 330%. 

The optimum pitch serves as the right trade-off between 
the attractive forces of the fins on the condensate and the 
channel area between the fins to drain the condensate from 
the tubs. A pitch smaller than the optimum has too narrow a 
gap between fins to efficiently allow condensate run-off. 
This is because the condensate drawn to adjacent fins only 
combines to create a thicker film as shown in Figure 5.8. 
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gure 5.8 Sketch of the effect of finning on condensation 
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VI. CONCLUSIONS RECOMMENDATIONS 

A, CONCLUSIONS 

1. The test apparatus successfully operates under 
vacuum conditions with no degrading effects caused by the 
presence of nnncondensable gases. 

2. The Sieder-Tate coefficient was deternined to be 

0.034 ± 0.001 using both an instrumented smooth tube and the 
modified Wilson plot technique. 

3. An oprimum fin pitch of 2.5 mm was found for the 
finned tubes tested. 

4. Heat-transfer enhancements up to 3007) were realited 
while the increase in the tuba outer surface area (due to 
finning) was only 88% over the smooth tube. 

B. RECOMMENDATIONS 

1. Install a second, more reliable pressure transducer 
to replace the manometer. 

2. Install a larger cooling-wat er sump to eliminate use 
of the once-through system from the fresh-water tap. 

3. Collect data using an insert for the test tube to 
decrease the internal convective resistance and, thereby, 
decrease the uncertainty in the calculated condensing coef- 
ficient. 

4. Take data for the finned tubes while varying fin 
height and thickness (as well as pitch). 
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APPMDIX A 

THERMOCOOPLE CALIBRATIOM 



1. EQUIPMENT USED 

a. Thermocouple Wire 

Copper-constant an, 0. 24 5-ram (0.01-in) Teflon- coated wire 
was used for all thermocouples. 

b. Calibration bath 

A Rosemont Engineering Model 91 3A calibration bath was 
used. A schematic representation of the barh is shewn in 
Figure A.1. 

1) Heating: Electrical 

2) Cooling: Liguid Nitrogen 

Note: Once a desired temperature is reached, the temp- 

erature is held constant by rapid cycling between heating 
and cooling. The bath is rated for teraperature fluctuations 
of less than 0.01®C. 

c. Thermocouple readout 

An HP 3054A Data Acquisition/Control System was used to 
obtain data. Resolution of the acquisition system was 1y<V. 

d. Hath temperature measurement 

A platinum resistance thermometer with an accuracy of 
0.01 C was used. 

2. PREPARATION 

a. Procedure 

The instrumented tube (with the wall thermocouples 
installed) and the steam thermocouples were immersed in the 
bath as well as the probes for the quartz thermometer. 

b. Analysis 

The computer program TCAL was used to monitor and store 
all thermocouple readings on a disk. A listing of the 
program is located in Appendix B. 
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3. CALIBRATION PROCEDORE 



a. The bath temperature was sat at about 10 ®C. 

b. When the bath temperature reached steady state, its 
value was entered into the computer. 

c. The computer automatically recorded and printed all 
thermocouple readings. 

d. The bath temperature was raised in increments of 10 
C to 90 °C and steps b and c were repeated for each incre- 
ment . 

4. CALIBRATION CURVES 

a. A least-squares method was used to generate a poly- 
nomial of the form: 

Dr = + a,.T + a^jT^ 

where: is the difference between the bath and thermo- 

couple temperatures, and 

T is the value of the thermocouples obtained using the 
seventh-order polynomial fitted for the Type T thermocouple 
wire used. This polynomial is listed in the beginning of 
the program TCAL. 

b. Coefficient values obtained were: 

a^ = 4.7338 X 10 -3 
a, = 7.6928 x 10*3 
a^ = -8. 0779 x 10 'S 

c. A plot of the curves is found in Figure A.1. The 
curve for the wall thermocouples gives a different reading 
due to the thermal conduction of the air temperature through 
the tube. The data used for calibration did not include 
this average value. Only the value obtained for the two 
steam thermocouples was used for all of the thermocouples. 
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Figure ft.1 Calibration curve for the thermocouples 
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APPENDIX B 
SAMPLE CALCOLATIONS 



1. HEAT LOSS 

A set of calculations was parforaed to estimate the heat 
loss of the apparatus dua to natural convection. These 
calculations were performed for two main sections; the 
boiler (with and without insulation» and the piping ro the 
test section. The minimum power input used during the thesis 
was 17 hW. 

a. Boiler (without insulation) 

1) Given tha following dimensions and properties of the 
toiler section: 



r = 0.152 4 m 
r^^jj ~ 0.1561 m 
L j = 0.4064 m 
k Q. =1.4 W/m-K 

2) Given the following temperatures; 

T;,53 = lOO^C 

T^.s= 95°C 
T.^ = 25 *C 

T^ = (95+25)/2 = 60 ^C 

3) Given the following properties of ai 
erat ure : 

= 28.74x10-3 W/m-K 



r at tha film temp- 
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Pr = 0.702 

f' - 19.21x10~‘ m2/s 

0 = 0.003 K'l 

4) Compute the Grashof Number for the boiler assuming 
laminar conditions (Ra < 10^): 

Grg = g(J(T, 8 “T^)L3/:r.2 

Gta = 3.7x108 

5) Compute the Rayleigh Number for rhe boiler: 

Rag = Gr Pr 

Rag = 2.6x108 (< 109) 

6) Determine the applicability of a flat-plate analysis 
according to the method of Sparrow and Gregg (Ref. 13] 

L /ro,6 < (0.15/TS)Gr»/9 

2.60 <7.36 

8) Compute the average Nusselt Number (via flat-plate anal- 
ysis) ; 

iTu = [0.825 + 0.387Raj i/‘/C 1 + (0.492/?r) 9/1* ]s/27 ]2 

■fu = 8 1.3 

9) Compute the external heat-transf er coefficient: 
ho = iiu k /L 

ho = 5.75 H/m2-K 

10) Compute the external convective thermal resistance: 

R = 1/(ho2Trro,„L^) 

R = 0.436 K/W 

11) Calculate the wall resistance: 
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I 




I 



I 



if 




Ru- = In (r^^o/r;,s) / 27 t)c^Lb 
R^ = 0, 007 K/W 

12) Neglecting the wall resistance, calculate nhe heat less 
from the boiler: 

Qg = )/(R *^w) 

Q 8 = 169.3 W 

t. Boiler (with insulation) 

Qg = 29.2 M 
c. Piping 

Q p = 65.0 W 

a. Total (worst case) 

Q = Q 3 ♦ Q p 

Q = 23U.3 W (<< 17000 W of input power) 
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APPENDIX C 

STEAH-CLBANINS PRDCEDDBE 

Note: The procedures listed here provide an excell enr 

method of cleaning the apparatus, but will cause the walls 
of the apparatus to heat to temperatures over 100 C (212 F) . 
Once this happens, several hours are required before the 
walls can cool down such that they won't superheat the steam 
generated during an experimental run. Therefore, use this 
method only prior to operating at atmospheric pressure or 
when the system is badly contaminated. 

1. Ensure the water level in the boiler is at leasx six 
inches above the upper ends of the heaters. 

2. Remove the thermopile from the inlet side of the 
test section. 

3. Energize the power supply and adjust the rheostat 
until the voltmeter reads about 100 V. 

4. Close the recirculating valve. 

5. Open the fill/drain valve. 

6. Once vapor begins exiting the fill/drain line, allow 
the sysxem to steam fcr several minutes. 

7. To maximize the steam flow through the test section, 
fully open the throttle to the auxiliary condenser. 

8. When done, simultaneously close the drain/fill valve 
and open the recirculating valve. 

9. Circulate cooling water through the test tube and 
check for the presence of any dropwise condensation. 



10. Reinstall the thermopile 
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APPBNDIX D 

SYSTEH STAHT-OP AND SHOT-DOHN PROCEDOHES 
To Start the system: 

1. Fill the boiler to at least six inches above the 
upper level of the heaters. 

2. Energize the air ejector for ten minutes. 

3. Circulate water through rhe cooling water sump by 
opening rhe inlet valve and activating the siphon. 

4. Energize the power to the boiler and adjust -he 
rheostat until the voltmeter reads 90 7 for vacuum opera- 
tion, or 170 7 for atmospheric pressures. 

5. Energize the circulating pumps and adjust the 
throttle for the desired flow through the tube. 

6. Open the valve to the auxiliary condenser and adjust 

the flow: about 15% for vacuum operation or 30% for atmos- 

pheric runs. 

7. Observe the steam temperature indicated on th efront 
of the data acquisition unit until the voltage is reached 
corresponding to the saturation temperature of the the 
desired operating pressure. 

8. Energize the air ejec-or for two minutes. 

9. Adjust the the flow rare through the auxiliary 
condenser to bring the apparatus to the desired operating 
pressure. 

To secure the system : 

1. Secure cooling water to the test tube. 

2. Secure power to the boiler. 

3. Open the vent valve. 

4. Secure the water supply to the auxiliary condenser. 
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APPENDIX E 

COMPUTER PROGRAM LISTINGS 

The following pages contain the computer program lisr- 



used for this 


t hesis: 


SIEDER 


(page 


61) 


WILSCN 


(page 


65) 


DR? 


(page 


71) 


TCAL 


(page 


80) 
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PROGRAM SIEDER 



1000! FILE NAME: SIEDER 

1010! DISK NUMBER: 12 
,1020! REVISED: November 23. 1983 

1030! 

1040 COM /Cc/ C<7) 

1050 DIM Emf (10) ,Tui<5) 

iOGO DATA 0. 10086091 .25727,94369.-787345.8295.78025595.81 
1070 DATA -9247486589. 6. 97688E1 1 , -2. G6192E13. 3. 94078E14 
1080 READ C(*) 

1090 Kcu-335 

1100 Di=.0127 

1110 Do=. 01905 
1120 Dr». 015875 
1130 L>. 13335 

1140 L1>. 060325 

1150 L2-. 034925 

1160 PRINTER IS 701 
1170 BEEP 
1180 CLEAR 709 

1190 INPUT "ENTER MONTH. DATE AND TIME <MM:DD:HH:MM:SS)".BS 
1200 OUTPUT 709:"TD”:BS 
1210 Series:! 

1220 OUTPUT 709:"TD" 

1230 ENTER 709: AS 

1240 PRINT USING "lOX. ""Month, date and time: "".14A":A$ 

1250 BEEP 

1260 INPUT "ENTER DISK NUMBER". Dn 
1270 PRINT 

1280 PRINT USING "1 OX .""NOTE : Program name: SIEDER 

1290 PRINT USING "16X.""Dislc number = "".DD":Dn 
1300 BEEP 

1310 INPUT "ENTER INPUT MODE ( I »305^A .2=F!LE)" , Im 
1320 IF Im»1 THEN 
1330 BEEP 

1340 INPUT "GIVE A NAME FOR THE DATA' FILE" . D_fi ieS 
1350 CREATE BDAT D_fileS.10 
1360 ELSE 
1370 BEEP 

1380 INPUT "GIVE THE NAME OF THE DATA FILE".D_f ileS 
1390 BEEP 

1400 INPUT "ENTER THE NUMBER OF RUNS STORED". Nrun 

1410 PRINT USING "16X.""This analysis was performed for file "" . 1 OA" : D_f i le$ 
1420 END IF 
1430 BEEP 

1440 INPUT "GIVE A NAME FOR PLOT DATA FILE". Plots 
1450 BEEP 

1460 INPUT "ENTER OPTION FOR END-FIN EFFECT ( 1 > Y . 0»N)" . Ife 

1470 IF lfe=0 THEN PRINT USING "16X,""This analysis neglects end-fin effect" 

1480 IF Ife=1 THEN PRINT USING "16X.""This analysis includes end-fin effect" 

1490 CREATE BDAT Plots. 5 

1500 ASSIGN SFile TO D_fileS 

1510 ASSIGN SFilep TO Plots 

1520 J=0 

1530 Sx=0 

1540 Sy-0 

1550 Sxs=0 

1560 Sxy>0 

1570 IF Im=1 THEM 

1580! READ DATA THROUGH THE DATA ACQUISITION SYSTEM 



1590! IF THE INPUT MODE (Im) = 1 
1600 BEEP 

1610 INPUT "ENTER FLOWMETER READING" ,Fm 
1620 OUTPUT 709:"AR AF20 AL30 VR1" 

1630 FOR 1=0 TO 10 
16A0 OUTPUT 709:"AS SA” 

1650 IF I>A THEN 

1660 Se=0 

1670 FOR K=0 TO 19 

1680 ENTER 709 :E 

1690 Se-Se+E 

1700 NEXT K 

1710 Emf<D=ABS(Se/20) 

1720 ELSE 

1730 ENTER 709;E 

1740 Emf(I)-ABS(E> 

1750 END IF 
1760 NEXT I 

1770 OUTPUT 713:"T1R2E" 

.1780 WAIT 2 
1790 ENTER 713;Ttt 
1800 OUTPUT 713;"T2R2E" 

1810 WAIT 2 

1820 ENTER 713;T2 

1830 OUTPUT 713;"T1R2E" 

1840 WAIT 2 
1850 ENTER 713:T12 
1860 T1=(T11+T12)-.5 

1870 ELSE 

1880! READ DATA FROM A USER-SPECIFIED FILE IF 

1890! INPUT MODE <Im) = 2 

1900 ENTER 9File:Emf(») ,T1 .T2.Fm 

1910 END IF 

1920 Tavg=<Tl+T2)*.5 

1930 Twall=0 

1940 FOR 1=5 TO 10 ^ 

1950 Tu<I-5)=FNTvsv(Emf(I)) 

1970 Twall=Twall+Tw<I-5) 

1980 NEXT I 

1990 Tuall=Twall/6 

2000 Cpw=FNCpijj( Tavg) 

2010 Rhow=FNRhow< Tavg) 

2020 Md=5.00049E-3+6.9861937E-3*Fm 

2030 Md*Md»(1 .0365-1 .96644E-3-T1 +5 .252E-6*TI '2)/. 995434 
2040 Mf=Md/Rhou 

2050 Vw-Mf/(PI«Dr2/4) 

2051 T2c=T2-( .0138+.001*Vui-'2) 

2060! T2c=T2-.004»Uw'2 

2070 0=Md»Cpu*(T2c-T1 ) 

2080 D tw=Q*LOG < Do /D i ) / ( 2*P I*Kcu»L ) * . 5 
2090 TuaIl = Twall-Dtu) 

2100 Lmtd»(T2c-Tt )/LOG((Tuall-Tt )/<TualI-T2c)) 

2110 Kw=FNKw< Tavg ) 

2120 P1=PI»<Di+Do) 

2130 P2«PI«<Di+Dr) 

2140 Al=TDo-Di)»PI«Do 
2150 A2=<Dr-Di)-PI»Dr 
2160 Hi=Q/(PI»Di»L*Lmtd) 

2170 IF If 0=0 THEN 
2180 Hic=Hi 
2190 GOTO 2300 
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2200 END IF 

22:0 M1=<Hi*P1/(Kcii*rt1 )) '.5 
2220 M2=(Hi*P2/<Kcu*A2))'-.5 

2230 Fel'FMTanhdIl-LI )/<H1*LI ) 

22A0 Fe2=FNTanh<M2*L2)/(M2*L2) 

2250 Hic=0/fPI*Di*(L+L1*Fei+L2»Fe2)*Lmtd) 

2260 IF ABS((Hi-Hic>/Hic)>.01 THEN 

2270 Hi=(Hic+Hi)*.5 

2280 GOTO 2210 

2290 END IF 

2300 PRINT 

2310 PRINT USING "1 OX , ""Posi tion number : I 2 3 4 5 

6 

2320 PRINT USING "1 OX . ••"Nal 1 temperature <Deq C) : .6(DD . DD , 1 X )” : Tu< *) 

2330? CALCULATE THE NUSSELT NUMBER 
2340 Nu=Hic«Di/Ku 
2350 Muu*FNMuw( Tavq ) 

2360 Re=Rhow*yw*D i/tluw 
2370 Cf = < Muw/FNMuii)< Twal 1 ) ) ' . 1 4 
'2380 Pru=FNPru< Tavg) 

2390 X=Re'.8*Prw*.3333*CF 

2400? COMPUTE COEFFICIENTS FOR THE LEAST-SOUARES-FIT 
2410? STRAIGHT LINE 
2420 OUTPUT 9Filep:X.Nu 

2430 PRINT USING "tOX .••"Tual I Tin Tout Lmtd Vw X Nu 

2440 PRINT USING "1 OX . 4<2D . 2D ,2X) , Z .DD ,2X . 50 . D . 2X . 4D . D" ; Tuial 1 . T1 . T2c . Lmtd . Vu . X . 
Nu 

2450 Sx=Sx>X 
2460 Sy=Sy+Nu 
2470 Sxs=Sxs+X»X 
2480 Sxy=Sxy+X»Nu 

2490? STORE RAW DATA IN A USER-SPECIFIED FILE IF 
2500? INPUT MODE dm) = 1 

2510 IF Im=1 THEN OUTPUT 9Fi le:Emf <*) , T1 . T2 .Fm 

2520 BEEP 

2530 J-J+1 

2540 IF Im=1 THEN 

2550 INPUT “ARE YOU TAKING MORE DATA ( 1 =YES . 0=N0 )?" ,Go_on 
2560 Nrun=J 

2570 IF Go_on=1 THEN 1570 
2580 ELSE 

2590 IF J<Nrun THEN 1570 
2600 END IF 
2610? Ci=Sxy/Sxs 

2620 Ci=<Nrun»Sxy-Sy*Sx)/<Nrun»Sxs-Sx '2) 

2630 Ac'*(Sy-Ci»Sx)/Nrun 
2640 PRINT 

2650 PRINT USING "I OX . ”"Sieder-Tate Coefficient - "" ,D . 4D" :Ci 
2660 PRINT 

2670 PRINT USING "1 OX . '•"Least-Squares Line: 

2680 PRINT USING "1 2X . '••’Slope = '••' ,MD .5DE . '' :Ci 

2690 PRINT USING ••12X.'”’Intercept = '••' ,MD . 5DE . : Ac 

2691 PRINT 

2700 IF Im=1 THEN 
2720 BEEP 

2730 PRINT USING '•10X,''”NOTE: '"', 11 ,"" data runs were stored in file ''•' ,8A” ;Nru 
r. . D_f i leS 

2731 ELSE 

2732 PRINT USING '•10X,''''NOTE: The above analysis was performed for file '''’.14A'' 
;D_f i leS 
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2740 

2750 

2760 

2770 

2780 

2730 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 



END IF 

PRINT USING “16X/*’Tlot data are stored in file -M 4A** ; P lots 
ASSIGN •iFile TO * 

ASSIGN ^Filep TO ^ 

END 

DEF FNRhou(T) 

Ro= 1 006, 35724-T*( ,774489- T^< 2. 262459E-2-T*3.03304E-4)> 

RETURN Ro 
FNEND 

DEF FNPru(T) 

Pnij=FNCpu< T ) ^FNHulu < T ) /FNKu( T ) 

RETURN Pru 
FNEND 

DEF FNriuu<T> 

A=247.8/<T+I33.15) 

Muu*2.4E-5*10''A 
RETURN Mow 
FNEND 

DEF FNKu(T) 

KU-.572183504477-H .52770 1 21 209E-3»T 

RETURN Ku 

FNEND 

DEF FNTvsv(Emf) 

COM /Cc/ C<7) 

S«iin*C ( 0 ) 

FOR 1*1 TO 7 
SujT)*Sujn+C ( I)*Emf 
NEXT I 
RETURN Stun 
FNEND 

DEF FNCpui<T) 

Cpui*<4,21 l20858-T*(2.26326E-3-T*<4.42361E-5+T*2.71428E-7)))*l000 

RETURN Cpu 

FNEND 

DEF FNTanh(X) 

P*EXP(X) 

0=EXP<-X) 

Tanh*(P+Q)/<P-Q) 

RETURN Tanh 
FNEND 
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PROGRAK WILSON 



1000! 


FILE MAHE: WILSCM 






1010! 


REVISED: December 5. 1983 






1020! 








1030 


COM /Cc/ C<7) 






1040 


DATA 0.10086091 .25727.94369.-767345.3295.78025595. 
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1050 


DATA -9247486589 . 6 . 97688E 1 1 . -2 . 661 92E 13 . 3 . 94078E 1 4 






1060 


READ C(») 






1070 


DIM Emf(4) 






1080 


130175 






1090 


LI >.060325 






1100 


L2-. 034925 






1110 


Do>. 01905 






1120 


Di=.0127 






1130 


Dr=. 015785 






1140 


Kca*335 






1 150 


Rm=Do'*L0G(Do/D i )/(2*Kca) 






1160 


PRINTER IS 701 






1170 


BEEP 






‘1180 


CLEAR 709 






1190 


INPUT "ENTER MONTH. DATE, AND TIME <MM:DD :HH :MM:SS" ,B$ 




1200 


OUTPUT 709;**TD”;BS 






1210 


Jp=0 






1220 


OUTPUT 709;”TD*' 






1230 


ENTER 709 :A$ 






1240 


PRINT USING "10X. ""Month, date and time : "".14A" 


:AS 




1250 


BEEP 






1260 


INPUT "ENTER DISK NUMBER", Dn 






1270 


PRINT 






1280 


PRINT USING "1 OX. ""NOTE: Program name : M WILSON" 






1290 


PRINT USING "16X,”"Disk number = "",DD":Dn 






1300 


BEEP 






1310 


INPUT "ENTER INPUT MODE ( 1 =3054A ,2=FILE)" . Im 






1320 


IF Im=1 THEN 






1330 


BEEP 






1340 


INPUT "GIVE A NAME FOR THE DATA FILE",D fileS 






1350 


CREATE BDAT D file$,10 






1360 


ELSE 






1370 


BEEP 






1380 


INPUT "GIVE THE NAME OF THE DATA FILE".D files 






1390 


PRINT USING "16X,""This analysis is for data in fi 


le "", 


, 14A":D_f ileS 


1400 


BEEP 






1410 


INPUT "ENTER THE NUMBER OF RUNS STORED". Nrun 






1420 


END IF 






1430 


BEEP 






1440 


INPUT "GIVE A NAME FOR PLOT-DATA FILE". Plots 






1450 


BEEP 






1460 


INPUT "ENTER OPTION ( 1 -OCT ,2=T-PILE . 3>AVE>" , I tm 






1470 


BEEP 






1480 


INPUT "ENTER OPTION FOR END-FIN EFFECT <1=Y,0=N)". 


Ife 




1490 


IF Itm»1 THEN PRINT USING ’M 6X /•"This analysis uses OCT 


readings 


1500 


IF Itm*2 THEN PRINT USING "16X,""This analysis uses T-PILE readings 


1510 


IF Itm’3 THEN PRINT USING "16X,""This analysis uses average of QCT and 


ILE readings 






1520 


IF Ife=1 THEN PRINT USING "l6X,""This analysis me 


lades 


end^fin effect 


1530 


IF lfe=0 THEN PRINT USING "16X,""This analysis neglects 


end-fin effect 


1540 


CREATE BDAT Plots, 10 






1550 


ASSIGN SFilep TO Plots 






1560! 


Cia*.040 






1570! 


Cil=.028 
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1580 


J j = 0 




1590 


Ci=.03 




1600 


J = 0 




1610? 


C i = < C i «JL+C 1 i > * . 5 




1620 


Sx = 0 




1630 


Sy = 0 




1640 


Sxs*0 




1650 


Sxy = 0 




1660 


PRINT 




1670 


PRINT USING OX /’-Iteration munber 




1680 


IF Jj=0 OR Jp=l THEN 




1690 


PRINT 




1700 


PRINT USING -12X/-T1 T2 Tsat 


Lmtd 


1710 


END IF 




1720 


ASSIGN ?Fiie TO D file$ 




1730 


IF Im-1 AND Jt=0 THEN 




1740? 


READ DATA THROUGH THE DATA ACQUISITION SYSTEM 


1750? 


IF THE INPUT MODE (Im) = 1 




1760 


BEEP 




1770 


INPUT -ENTER FLOWMETER READING”, Fm 




‘1780 


OUTPUT 709; -AR AF60 AL63” 




1790 


OUTPUT 709; -AS SA” 




1800 


EtD=0 




1810 


FOR 1=1 TO 20 




1820 


ENTER 709 ;Et 




1830 


Etp=E tp+Et 




1840 


NEXT I 




1350 


Etp=Etp/20 




1860 


OUTPUT 709; -AS SA” 




1870 


Ptran=0 




1880 


FOR 1=1 TO 50 




1890 


ENTER 709;Pt 




1900 


Ptran=Ptran+P t 




1910 


NEXT I 




1920 


Ptran=Pt ran/50 




1930 


OUTPUT 709; -AS SA” 




1940 


ENTER 709;Bvol 




1950 


OUTPUT 709; -AS SA” 




1960 


ENTER 709;Bamp 




1970 


OUTPUT 709;-AR AF20 AL24- 




1980 


FOR 1=0 TO 4 




1990 


OUTPUT 709; -AS SA” 




2000 


ENTER 709;Emf(I> 




2010 


Emf <I>=ABS<Emf<D) 




2020 


NEXT I 




2030 


OUTPUT 713;-T1R2E- 




2040 


WAIT 2 




2050 


ENTER 713;T11 




2060 


OUTPUT 713;-T2R2E- 




2070 


WAIT 2 




2080 


ENTER 713;T2 




2090 


OUTPUT 713;-T1R2E- 




2100 


WAIT 2 




2110 


ENTER 713;T12 




2120 


Ti=(T11+Tl2)*.5 




2130 


CLEAR 713 




2140 


ELSE 




2150? 


READ DATA FROM A USER-SPECIFIED FILE IF 


INPUT 


2160 


ENTER @Fiie;Bvoi ♦Bamp .Ptran,Etp.Emf(*),Fm,T1 ,T 


2170 


END IF 




2180 


Tsat“FNTvsv(( Eraf ^0)+Emf(l>)».5) 





■.DD”;Jj+l 






MODE <Im) 



Y 
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2190 Ti=FNTvsv< Emf (2 ) ) 

2200 Grad=FNGrad< ( T I +T2) • . 5) 

2210 To-Ti+ABS<Etp)/(l0*Grad)»1 .E+6 
2220 IF Jj=0 THEN 
2230 Er1=ABS(Ti-T1> 

2240 PRINTER IS 1 

2250 PRINT USING T1 = •*" .DD . 3D" : T1 

2260 PRINT USING Ti = "" ,DD ,DD" : Ti 

2270 IF Erl >.5 THEN 
2280 BEEP 

2290 PRINT "OCT AND TC DIFFER MORE THAN 0.5 C" 
2300 BEEP 

2310 INPUT "OK TO GO AHEAD ( I =Y . 0=N)7" .Olcl 



2320 END IF 

2330 PRINT USING DT <QCT) - "",Z.3D":T2-T1 

2340 PRINT USING DT (T-PILE) = "" .Z .3D" : To-Ti 



2350 IF Ok 1=0 AND Erl >.5 THEN 3600 
2360 Er2»ABS((T2-T1 )- ( To-Ti ) ) /( T2-T1 ) 

2370 IF Er2>.05 THEN 
2380 BEEP 

2390 PRINT "QCT AND T-PILE DIFFER MORE THAN 57." 

2400 BEEP 

2410 INPUT "OK TO GO AHEAD < 1=Y.O=N)?",Ok2 

2420 IF 0k2=0 AND Er2>.05 THEN 3600 

2430 END IF 

2440 PRINTER IS 701 

2450 END IF 

2460! CALCULATE THE LOG-MEAN-TEMPERATURE DIFFERENCE 

2470 IF Itm=1 THEN 

2430 Tf=T1 

2490 T1=T2 

2500 END IF 

2510 IF Itm=2 THEN 

2520 Tf=Ti 

2530 Tl-To 

2540 END IF 

2550 IF Itm»3 THEN 

2560 Tf=<T1+Ti)*.5 

2570 Tl=(T2+To)».5 

2530 END IF 

2590 Tavg-(Tf+Tl)-.5 

2600 Trise=Tl-Tf 

26 1 0 Lmtd'Tr ise/LOG( ( Tsat-Tf ) / ( Tsat-Tl ) > 

2620 Cpu=FNCpw<Tavg) 

2630 Rhou)=FNRhou< Tavg) 

2640 Ku=FNKw< Tavg) 

2650 Muuia=FNMuu( Tavg) 

2660 Prw=FNPrw( Tavg) 

2670 Mdt=5.00049E-3+6.9861937E-3*Fm 

2680 Md=Mdt*<1 .0365-Tf»(1 .96644E-3-Tf *5 .252E-6) )/ .995434 

2690 Vf=Md/Rhou 

2700 Vw=Vf/(PI-Di*2/4) 

2710 Trise = Trise- .004«Uw’'2 
2720 Q=Md»Cpw»Tr ise 
2730 Qp=0/(PI»Do*L) 

2740 Uo=Op/Lmtd 

2750 Re=Rhou*'/u»Di/Muua 

2760 Fe1=0 

2770 Fe2=0 

2780 Cf=1 

2790 Two=Tsat-5 
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2300 Tf ilm=Tsat/2+Two/2 
2810 Kf=FNK(j(Tf ilm) 

2820 Rhof “FNRhow( Tf i lm> 

2830 Haf»FNMuu<Tf ilm) 

2840 Hfgp=FNHfg<Tsat)+.68»FNCpw(Tf i Im) •< Tsat-Tuo > 

2850 New=Kf*(Rhof2«9.799»Hfgp/(l1af*Do*Qp))'' .3333 
2860 Ho=. 655-Mew 
2870 Tuoc=Tsat-Qp/Ho 

2880 IF ABS((Tuoc-Two)/Tuoc)>.001 THE.M 

2890 Two=Twoc 

2900 GOTO 2800 

2910 END IF 

2920 Cf=1.0 

2930 0n*ega=Re'.S-Prw'.3333-Cf 

2940 Hi=Kw/Di-Ci-Qmeqa 
2950 IF Ife-0 THEN 3040 
2960 P1=PI-<Di+Do) 

2970 P2=PI-<Di+Dr) 

2980 Al=<Do-Di)*PI-<Di+Do)-.5 
2990 A2-<Dr-Di>-PI*(Di+Dr)*.5 

3000 M1-<Hi-P1/(Kca-A1))'‘.5 

3010 t12-<Hi-P2/(Kca-A2))''.5 

3020 Fe1=FNTanh<M1-L1)/<l11*L1 ) 

3030 Fe2=FNTanh(M2-L2)/<M2-L2) 

3040 Dt-0/<PI»Di-(L+U-FeUL2-Fe2)-Hi> 

3050 Cf c*<t1uwayFNI1u.w< Tavg+Dt ) ) ' . 1 4 
3060 IF ABS(<Cfc-Cf)/Cfc)>.01 THEN 
3070 Cf-(Cf+Cfc)*.5 
3080 GOTO 2930 
3090 END IF 

3100 X=Do-New/<Omega*Kw) 

3110 Y=New-n/Uo-Rm) 

3120! COMPUTE COEFFICIENTS FOR THE LEAST-SOUARES-FIT STRAIGHT LINE 

3130 IF Jp=l THEN OUTPUT ^FileprX.Y 

3140 Sx=Sx+X 

3150 Sy=Sy+Y 

3160 Sxs=Sxs+X*X 

3170 Sxy-Sxy+X-Y 

3180? STORE RAW DATA IN A USER-SPECIFIED FILE IF INPUT MODE (Im) - 1 

3190 IF Im-1 AND Jj=0 THEN OUTPUT SFilerBvol .Ba/np.Ptran.Etp.Enf <•> ,Fm.T1 ,T2 

3200 IF Jj=0 OR Jp-1 THEN PRINT USING '*8X . 5<2X ,3D .DD) .2(2X . D .5D)" ; Tf , T1 . Tsat ,L p 

td.Uw.X.Y 

3210 BEEP 

3220 J-J+1 

3230 IF Im-1 AND Jj=0 THEN 

3240 INPUT "DO YOU HAVE MORE DATA ( 1 =Y . 0-N )?" . Go_on 
3250 Nran=J 

3260 IF Go_on-l THEN 1730 
3270 ELSE 

3280 IF J<Nru.n THEN 1730 
3290 END IF 

3300 Sl*<Nrun-Sxy-Sy-Sx)/<Nran-Sxs-Sx '2) 

3310 Ac-<Sy-Sl-Sx)/Nnin 

3320 Cic-l/Si 

3330 Jj*Jj+1 

3340 IF Jp=l THEN Jp=2 

3350 IF ABS<<Cic-Ci)/Cic)>.001 THEN 

3360 Ci=(Cic+Ci)-.5 

3370 PRINT USING "1 OX. ’"'Intermediate Sieder-Tate coefft = ’'”.Z.4D":Ci 
3330 GOTO 1600 
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tLL.^C 

3^00 IF Jp=0 THEN Jp-I 
3410 END IF 
3420 IF Jp-1 THEN 1 600 
3430 Ci=<Ci+Cic)».5 
3440 PRINT 

3450 PRINT USING ••10X/*"Sieder-Tate coefficient = •"•.Z.4D";Ci 

3460 PRINT 

3470 PRINT USING "1 OX .""Least-Squares Line: 

3480 PRINT USING "10X,"" Slope • . Z . 5DE . " ;S 1 

3490 PRINT USING "10X,"" Intercept = "" ,MZ . 5DE , " : Ac 

3500 PRINT 

3510 IF Iin=l THEN 

3520 BEEP 

3530 PRINT USING "10X.""N0TE: "", 11 ,"" data runs are stored in file "" 
f 1 leS 

3540 ELSE 

3550 PRINT USING " 1 0X . ""NOTE : Above analysis was performed for data in 
10A";D_f lies 
3560 END IF 

3570 PRINT USING "16X,""Plot data are stored in file "" , 1 OA" ; P lots 
3580 ASSIGN SFile TO • 

3590 ASSIGN SFileo TO * 

3600 END 

3610 DEF FNRhou(T) 

3620 Ro- 1 006 . 35724-T-l .774489- T«(2 .262459E-2-T-3 . 03304E-4) ) 

3630 RETURN Ro 
3640 FNEND 
3650 DEF FNPry(T) 

3660 Pr-j=FNCDu(T)»FNMuu<T)/FNKui(T) 

3670 RETURN Prw 
3680 FNEND 
3690 DEF FNMuu(T) 

3700 A=247.8/(T+I33. 15) 

3710 Mu=2.4E-5-10'‘A 

3720 RETURN flu 
3730 FNEND 
3740 DEF FNKu<T) 

3750 X=(T+273. 15)7273.15 

3760 Kui=-.92247+X*<2.8395-X*< 1 .8007-X*< .52577- . 07344*X) ) ) 

37-70 RETURN Ku 

3780 FNEND 

3790 DEF FNTvsv(Emf) 

3800 COM 7Cc/ C(7) 

3810 Sum-C(O) 

3820 FOR 1=1 TO 7 

3830 Sum=Sum+C< I)»Emf * I 

3840 NEXT I 

3850 RETURN Sum 

3860 FNEND 

3870 DEF FNCpu<T) 

3880 Cpw=(4.21120858-T»<2.26826E-3-T*<4.42361E-5+2.71423E-7)))«1000 

3890 RETURN Cpu 

3900 FNEND 

3910 DEF FNTanh(X) 

3920 P=EXP(X) 

3930 Q=EXP(-X) 

3940 Tanh=(P+0)/(P-Q) 

3950 RETURN Tanh 

3960 FNEND 

3970 DEF FNGrad<T) 



8A":J.D 
file "" 
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3980 

3990 

^000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 



COM /Cc/ C(7> 

Grad-37. 9853+. 104388*T 

RETURN Grad 

FNEND 

DEF FNFvst(T) 

F=Z|66.444+T*<7.09451-T-1 .65808E-2) 

RETURN F 

FNEND 

DEF FNHfg(T) 
Hfg=2477200-2450*(T-10> 

RETURN Hfg 
FNEND 
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PROGRAM DRP 



1000 ! 

1010 ! 

1020 ! 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

'1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1230 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1^60 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 



FILE NFME: DRP 

REVISED: November 18, 1933 



COM /Cc/ C(7) 

DIN EmfdO) 

D8TA 0. 10086091 .25727 , 94369 . -767345 . 3295 . 78025595 . 81 
DATA -9247486589. 6. 97688E+11 .-2.66192E^13.3.94078E^14 
READ C(*) 



Di=.0127 
Do=, 01905 
Dr=, 015875 
Dssp* . 1 524 



Inside diameter of test tube 

Outside diameter of test tube 

Outside diameter of the outlet end 

Inside diameter of stainless steel test section 



Ax=^PI"*Dssp '2/4-PI*Do*L 



L=. 130175 
L1-. 060325 
L2-. 034925 
Kcu=385 
Ci=.034 



Condensing length 
! Inlet end ’*fin length” 

! Outlet end **f in length” 

! Thermal conductivity of Copper 
! Sieder-Tate coefficient 
Rm=Do*L0G<Do/Di)/(2^Kcu> ! Wall resistance based on outside area 
PRINTER IS 701 
CLEAR 709 
BEEP 

INPUT ”ENTER MONTH. DATE AND TIME <MM:DD:HH:MM:SS)” .DateS 
OUTPUT 709:”TD”;DateS 
OUTPUT 709:”TD” 



ENTER 709;DateS 

PRINT ” Month, date and time :”;DateS 

PRINT 

PRINT USING ”10X.””N0TE: Program name : DRP 

BEEP 

INPUT ”EMTER DISK NUMBER”. Dn 

PRINT USING ”16X,””Dislc numoer = ””,DD”;Dn 

BEEP 

INPUT ”ENTER INPUT MODE ( 1 =3054A , 2-FILE) Im 

IF Im-1 THEN 

BEEP 

INPUT -GIVE A NAME FOR THE RAW DATA FILE” , D_f i leS 
CREATE BDAT D_fileS.15 
ASSIGN §File TO D_fileS 
BEEP 

INPUT ”ENTER GEOMETRY CODE ( I -FINNED .0-PLAIN)”.Ifg 
OUTPUT ^File:Ifg 
IF Ifg-0 THEN 
BEEP 

INPUT ”WALL TEMPERATURE MEASUREMENT ( 1 -Y , 0-N>” . lut 

ELSE 

BEEP 

INPUT ”ENTER FIN PITCH. WIDTH AND HEIGHT” .Fp . Fy , Fh 
END IF 

IF lfg=0 THEN OUTPUT «File;lLJt 
IF Ifg=1 THEN OUTPUT ^F i le :Fp .Fu ,Fh 
ELSE 
BEEP 

INPUT ”GIVE THE NAME OF THE EXISTING DATA FILE” . D_f : leS 
PRINT USING ”16X,””This analysis was performed for data in file 



MOA” 



7 1 







I 

I 

I 



r 1 iei 

1550 BEEP 

1560 INPUT "ENTER THE NUMBER OF RUNS STORED". Nrun 
1570 ASSIGN «File TO D_file$ 

1580 ENTER .SFilerlfg 

1590 IF lfg=0 THEN ENTER «File:Iut 

1600 IF Ifg=l THEN ENTER iFx le:Fp .Fw.Fh 

1610 END IF 

1620 IF lfg=0 THEN 

1630 BEEP 

1640 INPUT "WANT TO CREATE A FILE FOR Nr vs F < I = Y , 0=N) . Inf 

1650 ELSE 

1660 lnf=0 

1670 END IF 

1680 IF Inf=1 THEN 

1690 BEEP 

1700 INPUT "GIVE A NAME FOR Nr vs F FILE".MrfS 
1710 CREATE BDAT NrfS.2 
1720 ASSIGN ?Nrf TO NrfS 
•1730 END IF 
1740 BEEP 

1750 INPUT "ENTER OPTION ( 1=0CT.2=T-PILE.3=AVE>". Itn 
1760 BEEP 

1770 INPUT "ENTER OPTION FOR END-FIN EFFECT ( I =Y ,0=N)" . Ife 

1780 IF Itm=1 THEN PRINT USING "16X,""This analysis uses OCT readings 

1790 IF Itm=2 THEN PRINT USING "16X,""This analysis uses T-PILE readings 

1300 IF Itm*3 THEN PRINT USING "16X.""This analysis uses average of OCT and T-P 

ILE readings 

1810 IF Ife“I THEN PRINT USING "16X,""This analysis includes end-fin effect’ 

1820 IF lfe=*0 THEN PRINT USING ”16X,""This analysis neglects end-fin effect' 

1830 PRINT USING "16X,""Sieder-Tate coefficient = "",Z.4D":Ci 
1840 BEEP 

1850 INPUT "GIVE A NAME FOR PLOT DATA FILE" ,P_f i leS 

1860 CREATE BDAT P_file$.5 

1870 ASSIGN L»Filep TO P_fileS 

1880 IF Iut=1 THEN 

1890 BEEP 

1900 INPUT "GIVE A NAME FOR WALL TEMPERATURE FILE",WtfS 

1910 CREATE BDAT WtfS.5 

1920 ASSIGN yFilet TO WtfS 

1930 END IF 

1940 BEEP 

1950 INPUT "ENTER OUTPUT VERSION < 1=SH0RT.2=L0NG>" .lov 

1960! IF Im=1 THEN 

1970! OUTPUT SFileilfg 

1980! IF lfg=0 THEN OUTPUT «File:Iut 

1990! IF Ifg=I THEN OUTPUT SFile:Fp.Fw.Fh 

2000! ELSE 

2010! ENTER »File:Ifg 

2020! IF lfg=0 THEN ENTER «File:Iwt 

2030! IF Ifg=l THEN ENTER SF i le :Fp ,Fu . Fh 

2040! END IF 

2050 IF lfg=0 THEN 

2060 PRINT USING "16X,""Tube type : PLAIN 

2070 ELSE 

2080 PRINT USING "16X.""Tube type : FINNED 

2090 PRINT USING "16X.""Fin pitch, width, and height (mm): "" , DD . D ,2X . Z . DD .2X . Z 

.D0":Fp.Fw.Fh 

2100 END IF 

2110 J=0 

2120 IF lov- I THEN 
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I 



liliii 



2130 


PRINT 




2! 40 


IF Inf=I then 




2150 


PRINT USING ••10X/”*Data Vu Uo Ho 


2?60 


PRINT USING ••tOX/*” <m/s) (W/m' 


2-K)(W/m'2 


2170 


ELSE 




2t30 


PRINT USING ••10X/**‘Data Uu 


Uo 


2130 
• • 


PRINT USING "10X,"" t> (m/s) <W/m'2-K) 


2200 


END IF 




2210 


END IF 




2220 


Sx=0 




2230 


Sy'O 




2240 


Sxs=0 




2250 


Sxy = 0 




2260 


Repeat : ! 




2270 






2280 


IF Im-1 THEN 




2230 


BEEP 




‘2300 


INPUT ’’LIKE TO CHECK NG CONCENTRATION (I=Y.O 


2310 


BEEP 




2320 


INPUT "ENTER FLOWMETER READING" .Fm 


■2330 


OUTPUT 709:"AR AF60 AL63 VR5" 




2340 


OUTPUT 709:'*AS SA** 




2350 


ENTER 709:Etp 




2360 


OUTPUT 709: "AS SA" 




2370 


OtrairO 




2380 


FOR 1=1 TO 50 




2390 


ENTER 709 :Vt 




2400 


Utran=Vtran+Ot 




2410 


NEXT I 




2420 


Utran=Vtran/50 




2430 


OUTPUT 709: "AS SA" 




2440 


ENTER 709:Bvol 




2450 


OUTPUT 709: "AS SA" 




2460 


ENTER 709:Bamp 




2470 


IF lut=0 THEN OUTPUT 709:"AR AF20 


AL24 URr 


2480 


IF Iut=1 THEN OUTPUT 709;"AR AF20 


AL30 VRV 


2490 


IF lut=0 THEN Nn-4 




2500 


IF Iuit = 1 THEN Nn=10 




2510 


FOR 1=0 TO Nn 




2520 


OUTPUT 709: "AS SA" 




2530 


IF I>4 THEN 




2540 


Se=0 




2550 


FOR K=0 TO 10 




2560 


ENTER 709 :E 




2570 


Se=Se+E 




2580 


NEXT K 




2590 


Emf<I)=ABS(Se/10) 




2600 


ELSE 




2610 


ENTER 709 :E 




2620 


Emf (I)=ABS<E) 




2630 


END IF 




2640 


NEXT I 




2650 


OUTPUT 713:"T1R2E" 




2660 


WAIT 2 




2670 


ENTER 713:T11 




2680 


OUTPUT 713:"T2R2E" 




2690 


WAIT 2 




2700 


ENTER 713 :T2 





Qp Vv 

(W/m'2) (m/s>* 



Ho 

(W/m'2-K) 



Qp 

(H/m 2> 



\Mg 



F 

(m/s 
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I u uu I Pu r / 1 3 : ■* r 1 K2E'* 



2720 WAIT 2 
2730 ENTER 713:T12 
2740 T1=(T1UT12>*.5 

2750 IF Ng=0 THEN 2800 
2760 BEEP 

2770 INPUT ‘‘ENTER NANOMETER READING (HL »HR . HRW)’\H1 ,Hr .Hru 
2780 Phg=Hl+Hr 
2790 Pwater*Hr-Hru 
2800 ELSE 

2810 IF Ifg=l OR lut=0 THEN 

2820 ENTER ^Fiie;BvoKBaiTiD,Vtran.Etp.EfT>f (0) .Emf (1 ) .Emf (2) ,Emf <3) .Emf (4) .Fm . TK T 
2 . Phg t Puater 
2830 END IF 

2840 IF lfg=0 AND Iut=l THEN ENTER i^Fi le :Bvol »Bajnp .Vtran ,Etp ,Emf (» ) »Fn . Tl , T2 .Ph 
g *Puater 

2850 IF J=1 OR J*I0 OR J=20 OR J-Nran THEN 

2860 Ng=1 

2870 ELSE 

'2880 Ng=0 

2890 END IF 

2900 END IF 

2910 Tsteam*FNTvsv((Emf (0)+Enf(l))^.5) ? COMPUTE STEAM TEMPERATURE 

2920 Troom=FNTvsv<Emf (3) ) 

2930 IF Iut=1 THEN 
2940 Tum*0. 

2950 FOR 1*0 TO 5 

2960 Tu(I)*FNTvsv(Emf <1-^5) > 

2970 Tum*Tum+Tu( I ) 

2980 NEXT I 
2990 Tum=Tum/6 
3000 END IF 

3010 Tcon*FNTvsv(Emf ( 4 ) ) 

3020 Psat*FNPvst ( Ts teajn > 



3030 


Rohg* 1 3529- 1 22-* ( T room-26 . 85) /50 






3040 


Rouater=FNRhou( Troom) 






3050 


Ptest=< Phg»Rohg-Pyater*Rowater )-*S .799/1000 






3060 


Pmm*Ptest/133.322 






3070 


Pkm=Ptest»1 .E-3 






3080 


Pks=Psat*1 .E-3 






3090 


Pk t=FNPvsv( Vtran ) » 1 .E-3 






3100 


Tsat*FNTvsp(Ptest ) 






3110 


Vs t*FNVvst < Tsteajn) 






3120 


Ppng* (Ptest-Psat ) /P tes t 






3130 


Ppst= 1 -Ppng 






3140 


Mfng=l /n +18. 015/28. 97-Psat/(Ptest-Psat)) 






3150 


Vfng*Mfng/<1 .608- .608*Mf ng) 






3160 


Mf ng*Mf ng» I 00 






3170 


Vf ng=Vf ng»1 00 






3180 


BEEP 






3190 


IF Iov=2 THEN 






3200 


PRINT 






3210 


PRINT USING 'M0X/'“Data set nujnber 






3220 


PRINT 






3230 


END IF 






3240 


IF Iov*2 AND Ng=I THEN 






3250 

r VMM 


PRINT USING “10X/’“ P Psat 


Ptran Tmeas 


Tsat 


/• 

3260 


PRINT USING "lOX."" <nun) <lfPa> (IcPa) 


(kPa) <C) 


<C) 



Mass 



3270 



PRINT USING “lOX^SOD.DD.ZX) .2(3D.DD.2X) 



2 ( M3D . D . 2X ) ’* : Pmn ♦ Pkn , PVs . P^c t » Ts te 
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h 



an , Tsat , Vf ng ,Mf ng 

3280 PRINT 

3290 END IF 

3300 IF Mfng>.5 THEN 

3310 BEEP 

3320 PRINT 

3330 IF Im=1 THEN 

3340 BEEP 

3350 PRINT 

3360 PRINT USING '' 1 OX .'"'Energi ze the vacaum system 

3370 BEEP 

3380 INPUT "OK TO ACCEPT THIS RUN ( I =Y ,0=N)?" .Ok 
3390 IF 0k=0 THEN 
3400 BEEP 

3410 DISP "NOTE; THIS DATA SET WILL BE DISCARDED!? " 

3420 WAIT 5 

3430 GOTO 2280 

3440 END IF 

3450 END IF 

3460 END IF 

3470 IF Im=I THEN 

3480 IF Ifg=1 OR lut=0 THEN 

3490 OUTPUT ■^F i le ; Bvol , Bamp , Vt ran ,Etp , Emf ( 0 > ,Enf < 1 ) .Emf <2> ,Emf<3) ,Emf(4) .Fm.TI , 

T2.Phg.Puater 

3500 END IF 

3510 IF lfg=0 AND Iuit=1 THEN OUTPUT ?F i le: Bvo! , Bamp .Utran .Etp.Emf (*) .Fra . T1 , T2.P 
hq .Puater 
3520 END IF 

3530 IF lfg=0 AND Iu»t-1 THEN OUTPUT 5File!;Tu(*> 

3540? ANALYSIS BEGINS 
3550 Ti=FNTvsv(Emf(2)> 

3560 Grad-FNGrad<ai+T2)».5> 

3570 To=Ti+ABS(Etp)/<10-Grad)-I .E+6 
3580 Er1=ABS(Ti-T1 ) 

3590 PRINTER IS 1 

3600 PRINT USING T1 (OCT) = . DD .20*’ : T1 

3610 PRINT USING Ti <TC) = "" .DD .3D" ; Ti 

3620 IF Eri>.5 THEN 
3630 BEEP 

3640 PRINT "QCT AND TC DIFFER BY MORE THAN 0.5 C" 

3650 BEEP 

3660 INPUT "OK TO GO AHEAD ( 1 =Y .0-N)?" .Ok I 
3670 END IF 

3680 PRINT USING DT <GCT) = "".Z.3D":T2-T1 

3690 PRINT USIf-IG DT (T-PILE) = "" . Z .3D" ; To-Ti 

3700 IF 0kl=0 AND Ert>.5 THEN 5100 
3710 Er2=ABS((T2-T1 )-<To-Ti))/(T2-T1) 

3720 IF Er2>.05 THEN 
3730 BEEP 

3740 PRINT "QCT AND T-PILE DIFFER BY MORE THAN 57." 

3750 BEEP 

3760 INPUT "OK TO GO AHEAD < 1 =Y ,0=N)7" .0k2 

3770 IF 0k2=0 AND Er2>.05 THEN 5100 

3780 END IF 

3790 PRINTER IS 701 

3800 IF Itm=1 THEN 

3810 T1i=T1 

3820 T2o=T2 

3830 END IF 

3840 IF Itm=2 THEN 

3850 T1i=Ti 
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3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3340 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

41 10 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 



T2o=To 
END IF 

IF Itm=3 THEN 
II i-(T1+Ti)*.5 
T2o=<T2+To)*.5 
END IF 

Tavg= (T1 i+T2o)*.5 
Cpij=FNCou< Tavg) 

Rhow»FNRhow( Tavg) 

Md-5.00049E-3+6.9861937E-3*Fin 

Md*Md*( 1 .0365-1 .96644E-3^Tavg+5.252E-6»Tavg'2)/.995434 

Mf=Md/Rhou 

Vij=Hf/(PI»Di'2/4) 

T2o=T2o-< .0138+.001*^/u^2) 

Q=Md*Cpw*< T2o-T 1 i) 

0p*Q/(PI*Do»L) 

Kw=FNKw( Tavg) 

Muu=FNMuw( Tavg ) 

Rei =Rhow*V^u*D i /Muw 
Prw=FNPrw< Tavg) 

Fe1=0. 

Fe2*0. 

Cf = t. 

Hi “Ku#-»Ci/D i*Rei" .8*Prw ' . 33 33 *08 
Dt=Q/(PI*Di*(L+L1*Fe1+L2*Fe2)»Hi> 

Cfc=*(Miiw/FNMuu( Tavg+Dt ) ) ' . 1 4 

IF ABS((Cfc-Cf)/Cfc)>.01 THEN 

Cf=<Cf+Cfc)*.5 

GOTO 4090 

END IF 

IF lfe»0 THEN GOTO 4250 
Pl=PI»(Di+Do) 

A1*(Do-Di)*PI*<Di+Do)*.5 
M1=(Hi*P1/(Kcu*A1 )>" .5 
P2-PI*(Di+Dr) 

A2=(Dr-Di )»PI*(Di+Dr)* .5 
M2=(Hi*P2/(Kcu»A2)>" .5 
Fe1=FNTanh(M1»L1 )/<M1»L1 ) 

Fe2=FNTanh(M2*L2)/(M2»L2) 

Lmtd*(T2o-T1 i) /LQG( ( Tsteain-TI i ) / ( Tsteajn-T2o) > 

Uo»Q/ < Lm td*PI»Do*L) 

Ho=I/n/Uo-Do*L/(Di*<L+U*FeUL2*Fe2)*Hi)-Rm) 

Dtc=Q/(PI*Di»<L+L1*Fe1+L2*Fe2)*Hi) 

IF ABS((Dtc-Dt)/Dtc)>.01 THEN 4090 
Hf g=FNHf g( TsteajT)) 

Tuo=Tsteajr>-Qp/Ho 

Tf i lfn*Tsteajn/3+2*Tuo/3 

Kf = FNKui(Tf ilm) 

Rhof =FNRhou( T film) 
fluf *FNMuw( Tf i Im) 

Hpq*.651»Kf»<Rhof^2*9.81»Hfg/(Muf»Do*Qp)) ' .3333 

Y=Hpq*Qp".3333 

X*Qp 

Sx*Sx+X 

Sy=Sy+Y 

Sxs=Sxs+X"2 

Sxy=Sxy+X»Y 

OUTPUT $Filep:Qp.Ho 

Ql=500 ! TO BE MODIFIED 

Qloss=Ql/(100-25)*<Tsteam-Troom) ? TO BE MODIFIED 
Hf c=FNHf ( T con ) 



76 



4470 Hf-FNHf(T3team> 

4^80 Mdv=0 

4490? Bp=(Bvo 1-100) '2/5.7S ! BOILER POWER IM Watts 

4500 Bp = (8vol*100)''2/5.76 

4510 Mdvc=( (Bp-Qloss)-i1dv*(Hf-rifc) )/Hfo 

4520 IF ABS((Mdv-Mdvc)/Mdvc)>.01 THEN 

4530 Hdv=(Mdv+Mdvc)* . 5 

4540 GOTO 4510 

4550 END IF 

4560 Mdv= (Mdv+Mdvc) » . 5 

4570 Vg=FNVvst( Tsteajn) 

4580 Ov=Mdv*Vg/Ax 
4590 IF Inf“l THEN 

4600 F“(9. 799-»Do»Maf-»Hf g) /< Vv '2»Kf *< Tst earn- Two) ) 

4610 Nu=Ho-0o/Kf 
4620 Ret=Vv*Rhof »Do/Maf 
4630 Nr=Nu/Ret'.5 
4640 END IF 

4650 IF Inf=1 THEN OUTPUT SNrftF.Nr 
4660 IF Iov=2 THEN 

4670 PRINT USING "lOX."" T (Inlet) Delta-T 

4680 PRINT USING ”10X.”" QCT TC OCT T-PILE 

4690 PRINT USING ”1 OX .2<DD .DD ,2X> .2(Z .3D .2X)’* : T1 . Ti . T2-T1 . To-T i 
4700 PRINT USING •‘lOX,”” Vw Rei Hi Uo Ho q 

Vv 

4710 PRINT USING ”1 OX . Z . DD . 1 X ,5< MZ . 3DE . 1 X) . MZ . DD" ; Ow . Rei .Hi . Uo .Ho .Qp . Vv 

4720 END IF 

4730 IF Iov=l THEN 

4740 IF Inf=1 THEN 

4750 PRINT USING ”1 1 X .DD .2X . Z . DD .2X ,2<5D . D .2X) . Z .3DE . 1 X . Z . DD . 2< 1 X . 3D . DD ) " : J . Vw. 

Uo.Ho.Op.Vv.F.Nr 

4760 ELSE 

4770 PRINT USING "1 1 X .DD . 2X . Z . DD .2X .2<f1D . 4DE . 2X) , Z . 3DE .3X , Z . DD” : J . Vu . Uo . Ho .Op . V 

4730 END IF 
4790 END IF 
4800 IF Iro-1 THEN 
4810 BEEP 

4320 INPUT "WILL THERE BE ANOTHER RUN < 1=Y.0=N)?".Go_on 
4830 Nran»J 

4840 IF Go_on=1 THEN Repeat 
4850 ELSE 

4860 IF J<Nriin THEN Repeat 
4870 END IF 
4880 IF lfg=0 THEN 
4890 PRINT 

4900 S 1 = <Nrun*Sxy-Sy*Sx ) /(Nrun*Sxs-Sx'2) 

4910 Ac= (Sy-S 1 *Sx ) /Nran 

4920 PRINT USING "1 OX. ’’"Least-Squares Line for Hnu vs q curve; 

4930 PRINT USING "10X."" Slope = .MD .4DE" :S1 

4940 PRINT USING "10X,”" Intercept » "" ,MD . 4DE" : Ac 

4950 END IF 

4960 IF Im=1 THEN 

4970 BEEP 

4980 PRINT 

4990 PRINT USING "10X.”"NOTE: '”’. 11 ,"" data runs were stored in file "".10A":J. 

D_f ileS 

5000 END IF 

5010 BEEP 

5020 PRINT 

5030 PRINT USING ” 1 0X , '"’NOTE : ””.ZZ.”” X-Y pairs were stored in plot data file 
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ileS 

50^0 IF Inf=! THEN 

5050 PRINT USING ”16X,ZZ,'*“ pairs of Nr-F are stored in file ”** , 1 4A’* : J .Nrf S 

5060 END IF 

5070 ASSIGN ?Fiie TO * 

5030 ASSIGN «Fiiel TO ♦ 

5090 ASSIGN ^Fiiep TO * 

5100 END 

5110 DEF FNPvst < Tsteam) 

5120 DIM K(8) 

5130 DATA -7.691 234564 , -26 . 08023696 .-168.1 706546 , 64 . 23235504 .-118. 9646225 
5140 DATA 4.16711732.20.9750676 J .£9.6 
5150 READ K(^) 

5160 T=<Tstean+273. 15)7647.3 

5170 Sum^O 

5180 FOR N=0 TO 4 

5190 Sam«Su/n+K(N>*< 1-T)"<N+1 ) 

5200 NEXT N 

5210 Br=Sam/<T*( 1+K(5)*( 1-T)+K<6)^< 1-T) '2) )-< 1-T) /(K< 7)-n -T) *'2+K(3) ) 

5220 Pr=EXP<Br) 

5230 P=22120000*Pr 

5240 RETURN P 
5250 FNEND 
5260 DEF FNHfg(T) 

5270 Hfg=2477200-2450*<T-10) 

5280 RETURN Hfg 
5290 FNEND 
5300 DEF FNMuu(T) 

5310 A-247.3/<T+133.15) 

5320 Mu=2.4E-5*10''A 

5330 RETURN Ma 
5340 FNEND 
5350 DEF FNVvst<Tt> 

5360 P=FNPvst(Tt) 

5370 T*Tt+273.15 

5380 X=1500/T 

5390 F1=1/< l+T*1 .E-4) 

5400 F2-(1-EXP<-X))'^2.5*EXP(X)/X" ,5 

5410 B- . 00 1 5*F 1 - . 000942*F2- . 0004882*X 
5420 K=2»P/(461 .52*T) 

5430 V»<1+<1+2^B*K>".5)/K 

5440 RETURN V 
5450 FNEND 
5460 DEF FNCpu(T) 

5470 Cpui=4.21 120858-T*(2.26826E-3-T*(4.42361E-5+2.71428E-7*T)> 

5480 RETURN Cpu*l000 

5490 FNEND 

5500 DEF FNRhou(T) 

5510 Ro=999.52946+T»< . 0 1 269-T*< 5 . 4825 1 3E-3-T* 1 .234147E-5)) 

5520 RETURN Ro 
5530 FNEND 
5540 DEF FNPrw<T) 

5550! Pru»10'< 1 . 09976605-T*< 1 . 37^9326E-2-T*< 3 . 968875E-5-3 . 45026E-7*T) ) ) 

5560 Pru*FNCpu< T)*FNMau< T)/FNKw< T ) 

5570 RETURN Pru) 

5580 FNEND 
5590 DEF FNKu<T) 

5600? Ku=.5625894+T*<2.2964546E-3-T»< 1 . 509766E-5-4 . 0581 652E-8*T) ) 

5610 X=<T+273. 15)7273. 15 

5620 Ku»-.92247+X*(2.8395-X*< 1 .8007-X*< . 52577- . 07344*X) ) ) 

5630 RETURN Ku 



78 




J 

i 



I 

I 



5640 FMEND 

5650 DEF FNTanh(X) 

5660 P-EXP(X) 

5670 O-EXP(-X) 

5680 Tanh’(P+Q)/<P-Q) 

5690 RETURN Tanh 

5700 FNEMD 

5710 DEF FNTvsv<V) 

5720 COM /Cc/ C(7) 

5730 SuLm*C(0) 

5740 FOR 1=1 TO 7 
5750 Siim*SuLm+C( I 

5760 NEXT I 

5770? T*U»< .02635206856*-y^<9.735t3t3E-7-V-6.576805E-n )) 

5780 RETURN Sum 

5790 FMEMD 

5800 DEF FNHf(T) 

5810 Hf=T^<4.203849-T^(5,88132E-4-T*4.55l603t 7E-6)) 

5820 RETURN Hf^lOOO 

S830 FNEND 

5840 DEF FNGrad(T) 

5850 Grad=37.9853+. t04388^T 

5860 RETURN Grad 
5870 FNEMD 
5880 DEF FNTv3p<P) 

5890 Tu=nO 
5900 Ti-10 

5910 Ta=<Ta+Tl)-^.5 

5920 Pc=FNPvst<Ta) 

5930 IF ABS(<P-Pc)/P)>.00t THEN 

5940 IF Pc<P THEN Ti=Ta 

5950 IF Pc>P THEN Tu=Ta 

5960 GOTO 5910 

5970 END IF 

5980 RETURN Ta 

5990 FNEND 

6000 DEF FNPvsv(U) 

Goto P=8I 33. 5 133+2.23605 I E+4*V 
6020 RETURN P 
6030 FNEND 
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PROGRAM TCAL 



'00 ! ^ILE OfiME: TCAL 
110 ! REVISED: Decemoer 11. 1983 

! 20 ! 

130 ' COM /C;/ C(7> 

MO DIM Emf<10) .1(10) .DOO> 

150 DATA 0.10086091.25727,94369,-7573^5.3295.73025595.81 
160 DATA -9247486589, 6. 97688E1 1 .-2. 66192E13. 3. 94078EM 
170 READ C(*) 

180 CLEAR 709 

190 BEEP 

200 INPUT "ENTER MONTH, DATE AND TIME (MM :DD :HH :MM:SS>" , BS 
210 J=0 

220 OUTPUT 709:"TD":BS 

230 OUTPUT 709:"TD" 

240 ENTER 709: AS 

250 PRINT USING " 1 OX , ""Month . date and time = "’■,14A";AS 
260 BEEP 

270 INPUT "ENTER INPUT MODE (M3054A, 2=FILE)".Im 

280 IF Im=1 THEN 

290 BEEP 

300 INPUT "GIVE A NAME FOR DATA FILE" , D_f i leS 
310 CREATE BOAT D_fileS.5 
320 ELSE 

330 BEEP 

340 INPUT "GIVE NAME OF EXISTING FILE" , D_f i leS 
350 BEEP 

360 INPUT "ENTER NUMBER OF DATA RUNS STORED", Nran 
370 END IF 

380 BEEP 

390 INPUT "GIVE A NAME FOR PLOT FILE" ,P_f i leS 

400 CREATE BDAT P_fileS.5 

410 ASSIGN iPlot TO P_fileS 

420 ASSIGN ^File TO D_fileS 

430 IF I.m-1 THEM 

440 BEEP 

450 INPUT "ENTER BATH TEMPERATURE" , T_bath 
460 OUTPUT 709;"AR AF20 AL30" 

470 FOR 1=0 TO 10 

480 OUTPUT 709; "AS SA" 

490 ENTER 709:Emf(I) 

500 NEXT I 

510 OUTPUT 713:"T1R2E" 

520 WAIT 2 

530 ENTER 713:T1 

540 OUTPUT 713;"T2R2E" 

550 WAIT 2 

560 ENTER 713 :T2 

570 OUTPUT fFi le : T_bath . tmt ( » ) , T 1 , T2 
580 ELSE 

590 ENTER SFiie:T_bath,Emf (*) ,T1 ,T2 

600 END IF 

blO J=J+I 

620 Dua=0 

630 FOR 1=0 TO 10 

640 T(I)=FNTvsv(ABS(Emf (I) ) ) 

650 D( I) = T_bath-T< I ) 

660 IF I>4 THEN Dua=Dua+D(D 
670 NEXT I 

630 Dua=Dua/G 
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690 Dsa= < D < 0 ) +0( ! ) ) ♦ . 5 

700 OUTPUT lot : T_bath . Dsa, Dua 

710 PRINT 

720 PRINT USING ’MOX/^'Data set nuinber = 

730 PRINT USING •M0X/**’Bath T (C> QCT«1 cC) QCT-2 (C) 

740 PRINT USING OX . 3< 3D . 3D . 7X)** ; T_bath . T i , T2 

750 PRINT USING ”1 OX . "’’ThermocoiiP le readinas <Deq C) :'*’*** 

760 PRINT USING ’M OX ,6< 3D .DD ,3X) J6X'* ; T( -> 

770 PRINT USING 1 OX . ’'‘‘D iscrepanc les (Deg C) : 

730 PRINT USING H X . 6< (1Z . DD . 4X ) J 5X’’ ; D( - ) 

790 BEEP 

300 IF Im=1 THEN 

810 INPUT “ARE YOU TAKING HORE DATA < 1 - Y . 0 = N>?*\Go_on 
320 IF Go_on-1 THEN 430 
830 ELSE 

840 IF J<Nrun THEN 430 
850 END IF 
860 PRINT 
870 IF Im=1 THEM 

380 PRINT USING 'M OX . *'**N0TE : '*'*,DD/*** data sets are stored in file '*** . 1 4A'* : J .D 
1 leS 

890 ELSE 

900 PRINT USING ’* 1 OX , ’*'*NOTE : Above analysis uas oer’^ormed from file **** . 1 4A** : D_ 
files 

910 END IF 

920 PRINT USING *' 1 6X /’“Plot data are stored in file *”/ 1 0A“ : P_f i leS 
930 ASSIGN '?File TO * 

940 ASSIGN iPlot TO 
950 END 

960 DEF FNTvsv<Emf) 

970 COM /Cc/ C(7) 

980 Sam*C(0) 

990 FOR 1=1 TO 7 
1000 Siim=S»jjn+C( I ) ' I 
1010 NEXT I 
1020 RETURN San 
1030 FNEND 
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